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7.1 Introduction 

This chapter deals with X-ray emission from isolated neutron stars for which 
the energy for the observed X-rays is thought to originate from the rotation of the 
neutron star, or from an internal heat reservoir following formation. Rotation power 
, can manifest itself as pulsed emission, or as nebular radiation produced by a rela- 

' tivistic wind of particles emitted by the neutron star. Residual heat of formation 

is observed as soft X-ray emission from young neutron stars. Such thermal radia- 
I tion, however, can also be produced as a result of reheating from internal or external 

' sources. Rotation-powered pulsed and nebular X-ray emission, as well as thermal 

(«~^ I emission, can often be observed in a single object simultaneously; this is both fasci- 

nating and annoying, as one invariably contaminates the study of the other. There 
are also a handful of neutron stars for which the origin of the observed X-ray emission 
is unclear but may be related to the above processes; we will discuss those as well. 
2 ' Rotation-powered neutron stars are generally referred to as "radio pulsars" since 

^1 it is at radio wavelengths that the vast majority of the catalogued* population (cur- 

■ rently numbering ~1400) is observed. However, the radio emission is energetically 

unimportant, and we now know of several rotation-powered neutron stars that are 
not detected as radio so urces in spite of deep searches (e.g. Crawford et all 19981: 



I\jj . iMcLaughlin et al. I l2nmh . We therefore use the more physically motivated term 



' "rotation-powered." 

The total available spin luminosity in a rotation-powered pulsar is given by the 
rate of loss of rotational kinetic energy, E = Iujlo = Att^IP/P^, where / is the stellar 
moment of inertia and uj = 27r/P is the angular frequency with P the spin period. 
Thus a simple measurement of P and P for an isolated neutron star determines 
the available rotational power in a model-independent way, assuming a value for /, 
typically taken to be 10**^ g cm^. Also generally inferred for these sources are their 
surface di-polar magnetic fields, B = 3.2 x lO^^(PP)^/^ G (for B on the equator), 
and their characteristic spin-down ages, via Tc — P/2P. Although these latter two 
inferences are model-dependent {B because it assumes simple magnetic dipole brak- 
ing in a vacuum - which is almost certainly not the case - and Tc because it represents 
the true physical age only for the same assumption and also for a negligible initial 



* http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 7.1. P-P diagram for the 1403 currently catalogued rotation-powered pulsars. 
The 66 X-ray detected sources are indicated with an "X." Solid lines show constant 
E, and dotted lines show constant inferred surface di-polar magnetic field. "X- 
ray-detected" means pulsed, unpulsed or nebular emission, be it thermal or non- 
thermal. X-ray detected sources are generally those with the greatest spin-down 
luminosity, although the correlation is not perfect because of the wide range of 
source distances and observational selection effects. The indicated sources are 
summarized in Tables mi 17.21 and 17.51 Note that the X-ray-detected millisecond 
pulsars in the globular cluster 47 Tucanae do not have observable intrinsic P's 
due to contamination by acceleration in the cluster potential; for those sources 
(num bering 15), we have used inferred P's (see lFreire et all 120011 : IGrindlav et al.L 



spin period), the important point is that P and P provide, in addition to E, at least 
estimates of other important physical information. 

It is for this reason that P-P diagrams are so useful. Figure ITtI shows the P-P 
phase space, with all catalogued rotation-powered neutron stars with measured Ps 
indicated as dots. Lines of constant E and B are indicated. Note the two different 
populations of rotation-powered neutron stars: those with periods between ~16 ms 
and ~8 s, having B > 10"'^^ G, and the lower magnetic field {B < 10^° G), very 
rapidly rotating class {P ^ 100 ms), the millisecond pulsars. The division is made 
particularly obvious when binarity is considered, as the latter sources are mostly in 
binaries, while the former are only rarely in binaries. The millisecond pulsars are 
often referred to as "recycled" pulsars as they are thought to be formed in a past 
episode of mass transfer as part of an accreting X-ray binary phase (see chapter 1 
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by Psaltis). Finally, note the absence of sources having B > 10^'' G. Neutron stars 
having fields higher than this are thought to be "magnetars," discussed in chapter 
14 by Woods et al. 

All sources which have been detected in X-rays are identified with an "X" in 
Figure 17.11 For the purposes of this Figure, any X-ray detection, be it pulsed, 
unpulsed or nebular, thermal or non-thermal, is counted. Note the clear correlation 
with E ~ most X-ray identified rotation-powered neutron stars have large E and vice 
versa. The correlation is not perfect; the main reason for this is that the distances 
to these sources vary by over two orders of magnitude. Distant high-_B sources are 
harder to detect, and nearby \ow-E sources may be detected through their thermal 
emission. Note also the large number of millisecond pulsar detections is due to many 
of them being in X-ray studied globular clusters having larger scale heights, hence 
smaller line-of-sight absorbing column densities. 

Tables 17.11 17.21 and 17.51 provide a summary list of all X-ray-detected rotation- 
powered neutron stars, along with their most important properties. Note that the 
number of X-ray detected rotation-powered pulsars has more than doubled in the 
past 5 years (from 27 in the 1997 review by Becker & Triimper, to 66 today). This 
is due to improvements in X-ray sensitivity, in particular from the Rossi X-ray Tim- 
ing Explorer (RXTE), the Chandra X-ray Observatory and the XMM-Newton Tele- 
scope, and also improvements in radio telescope sensitivity, most notably at the 
Parkes, Arecibo, and Green Bank radio telescopes, which have permitted discoveries 
of sources that could later be followed up with X-ray observations. 



7.2 Magnetospheric Emission 

7. 2. 1 Motivation 

After the discovery i Staelin fc Reifenste in'.'l968^'Comella et al.Vl969') of the 



Crab pulsar in 1968 and the realization llOstriker fc Gunn. . .l969»:.Finzi fc Woli . .l96d) 
that the spin-down energy loss of the neutron star, as estimated from its measured 
P and F, could power all the visible radiation of the Crab nebula, it was evident 
that neutron stars are energetically important sources of particle production and 
acceleration. It is now believed that the radiating particles in pulsar wind nebulae 
(see H7.6|) are electron-positron pairs produced by electromagnetic cascades within 
the pulsar magnetosphere. Since the pulsar magnetosphere is also believed to be 
the origin of the pulsed non-thermal emission, understanding this emission will lead 
us to an understanding of the energetics of these sources and ultimately how the 
pulsar converts its spin-down power into both the pulsed and unpulsed observable 
radiation. However, the pulsed radio emission, which first enabled the discovery of 
pulsars, makes up a tiny (~ 10^^) portion of the energy budget. Pulsed emission 
at high energies, although harder to detect, is energetically much more significant 
(making up ~ 10"'^ of the spin-down power in the Crab and up to 0.1 in other pulsars) 
and thus is a more direct probe of the particle acceleration in the magnetosphere. 

Aside from the important issue of energetics, we also want to study the non- 
thermal emission at high energies in order to understand the physical processes 
in the pair cascades which generate the electron-positron pairs. It is widely be- 
lieved that generation of a pair plasma is required for the instabilities that produce 
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coherent radio emission in pulsars I MelroseL 200Cl|) . Why radio pulsars die after 
they have spun down to periods of several seconds may possibly be understood as 
their inability to produce pair plasmas when particle acceleration becomes too feebl e 



tnen mabinty to proauce pan plasmas wnen particle acceleration becomes too leebl t 
I Ruderman fc SutherlandL Il975t Iffibschman fc AronsL l200lt iHarding et all l2002a|) 



Studying the radiation and pair cascades at high energy is also a study in the phys- 
ical processes, such as one-photon pair creation and photon splitting, that occur 
only in the strong magnetic fields of neutron stars. Finally, the generation of non- 
thermal emission is intimately linked to the polar-cap heating which contributes to 
the thermal emission, particularly in older pulsars. 



7.2.2 Summary of Observations 

At X-ray energies, 66 rotation-powered pulsars have been detected, nearly 
half of whic h have periods P < 6 ms . Many of th ese detection s were first made 

but the number 



by ROSAT ifRecker fc Triimnerl l1997^ and ASCA 



has increased r apidly in recent year s with new observations by RXTE, Chandra and 
XMM-Newton l|Becker fc Aschenbac h. 2002). Table 1.1 hsts the pulsars thought to 
show non-thermal emission components in their X-ray spectra with observed pa- 
rameters, including detected emission at other wavelengths. The vast majority are 
also radio pulsars; many of the target sources were chosen from radio pulsar cat- 
alogs. Howe ver, there have been some detections of X-ray pulsars, such a s PSR 
J0537-6910 ("MarshaU et al.', 'l998f), PSR J1846-0258 ("Cotthclf et all 120 00*). PSR 
J1811-1926 (,Ib rii_cl al,, IMP , and PSR J0205-f6449 (Murra y et all 12002.) . with- 
out a radio epheme ris. Radio pulses hav e since been detected from only one of these 
(PSR .10205+6449: ICa,mi1o et all [ 2002(1 . Clearly, short period pulsars seem to have 
the highest levels of non-thermal X-ray emission. Most of the field millisecond pul- 
sars have purely non-thermal spectra. On the other hand, the millisecond radio 
pulsars in th e globular cluster 47 T ucanae that have been detected as point sources 



by Chandra (|Grindlav et all l2002() have spectra that are quite soft, and most likely 
thermal. The longer-period pulsars with detected X-ray emission have spectra that 
are dominated by thermal emission. 

7.2.2.1 Spectra 

Generally, the X-ray spectra of spin-powered pulsars show a mix of ther- 
mal and non-thermal components. Often several thermal components (discussed 
in M7.3|l and sometimes even several non-thermal components can be identified. 
Cleanly separating the non-thermal, or power-law, components from the thermal 
components can be difficult, especially in sources where the thermal components 
dominate (e.g. Gem inga, PSRs 0656-fl4 and B1055-52: fHalpern fc Wall Il997t 
Greiveldinger et all [l996). Detectors with sensitivity at energies up to at least 10 
keV, such as RXTE, Chandra and XMM, have been able to make the best measure- 
ments of non-thermal emission components. This non-thermal emission probably 
originates from the radiation of particles accelerated in the pulsar magnetosphere. 
Models for this emission will be discussed below. 

All of the 7-ray pulsars that have been detected by the Compton Gamma-Ray 
Observatory (CGRO) are also X-ray pulsars with measured non-thermal emission 
components. One interesting question is whether the non-thermal X-ray spectrum 
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Table 7.1. Non-Thermal X-ray Detected Rotation-Powered Pulsars 



is part of the same emission as that detected at 7-ray energies. For the youngest 
of these pulsars (Crab, PSR B1509— 58), the total power in pulsed emission peaks 
in the hard X-ray band (Figure I7.2|l . There is strong emission through the entire 
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Fig. 7.2. Broadband pulsed spectrum l|KuiDer et all 1200 If) of the Crab, as an 
example of emission characteristics of a young pulsar. 



X-ray and 7-ray bands and a smooth connection between the two. However, in the 
case of the Crab the s moothness of the sp ectrum is misleading as the situation may 
be more complicated. iKuiper et al. I ll200]J) have argued for several separate emission 
components in the Crab optical to 7-ray spectrum, evidenced by the strong frequency 
dependence of the interpeak emission and the peak 2 to peak 1 ratio, both of which 
have a maximum around 1 MeV. The middle-aged pulsars have comparatively weak 
non-thermal emission in the X-ray band since their power peaks at GeV energies, 
and there is a gap in the detected spectrum between the X-ray and 7-ray bands 
f Figure [7^ . In several cases (Vela, PSR B1055— 52) an ex trapolation bet \ Yeen t he 
two is plausible, but in others (Geminga, PSR B1706-44; ICotthelf et all 120021) a 
connection is not clear. In the case of Vela f Fieure 17.311. there is str ong evidence for 
two separate non-thermal X-ray components l)Harding et al.Ll2002b|) . one connecting 
to the 7-ray (OSSE) spectrum and the other possibly connected with the optical 
spectrum. 

Millisecond pulsars seem to be bright X-ray emitters, as a relatively large number 
of the known millisecond radio pulsars are X-ray pul sars. Most have hard power-law 
spectra with photon indices 1.5-2.0 (,Becker fc Aschenbach. .2Q0 ; ^'l , with the exception 



spectra with photon indices 1.5-2.0 (,Becker fc Aschenbach. .2Q0 ; ^'l , with the exception 
of a few (e.g. PSRs J0437-4715, .12124-8358: IZavlin et all l2002t ISa,kiira,i et all 



20011) which also have thermal components. Since they are too old to have detectable 



cooling emission, the presence of pulsed thermal emission from millisecond pulsars 

probably requires polar-cap heating. 

There are a number of n ewly discovered X-ray pulsars (PSR J1420— 6048jRoberts et al 



2001 



2001 



(PSR J1930-hl852; 
), (PSR J1105-6107; 



Camilo_etal 

Kaspi et all 



J2002a'), (PSR ,12229-^6114: iHalpern et alT 



1997) that lie in or near the error circles of 



EGRET unidentified T-rav sources |Hartman et al.lll999() . Since these pulsars were 
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Fig. 7.3. Broadband pulsed spectrum ijHarding et al.L l2002bfl of the Vela pulsar, 
as an example of emission characteristics of a middle-aged pulsar. 



discovered after the end of the CGRO mission and EGRET detected too few photons 
to search for pulsations, it is not known whether they are also 7-ray pulsars. There 
are also some X-ray pulsars discovered earlier (PSR B1823— 13, PSR B 1800— 21 



that are associated with EGRET sources and are candidate 7-ray pulsars (jKanbacl 
2OO2I) . 



7.2.2.2 Pulse Profiles 

The pulse profiles of the non-thermal X-ray emission usually show much 
more modulation and narrower pulses than does the thermal emission. The pulsed 
fraction for many pulsars with mixed emission rises from around 10% - 20% in the 
(0.1 - 1 kcV) band where the therm al emission dominates , to 80 % - 90% where 
the non-thermal emission dominates I Becker &: AschenbacbL 2002() . There is often 



an accompanying phase shift between the thermal and non-thermal pulses of the 
middle-aged pulsars. 

There seems to be no pattern in the relative phases of the non-thermal X-ray, 
radio and gamma-ray pulses. The non-thermal X-ray and radio pulses of many 
pulsars that are young or have very short periods, including the Crab (with the 
exception of the radio pre-cursor), PSRs J1617-5055, J0205+6114, and J1930-)-1852 
and most of the millisecond pulsars, are in phase and/or the profiles look very similar. 
Somehow for these pulsars, the radio and high-energy emission must originate from 
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Fig. 7.4. Multiwavelength pulse profiles fKanb achl l2002fl of the Crab, as an exam- 
ple of emission characteristics of a young pulsar. 



the same region in the magnetosphere. However, the profiles of these pulsars vary 
from double-peaked to both broad and narrow single-peaked shapes, and a few young 
X-ray pulsars (PSRs J1811-1926, J0537-6910 and J1846-0258) have no detected 
radio counterparts. The non-thermal X-ray pulses of middle-aged pulsars, on the 
other hand, are not in phase with their radio pulses, indicating a different location 
in the magnetosphere of the two emission components (e.g. Figure ^3}. The same 
trend seems to be true for the X-ray and 7-ray pulses of the 7-ray pulsars, with 
the p rofiles of the fastest, including millisecond (PSR J0218-I-4232, iKiiip er et al.1 
2000(1 ■ pulsars being alike and in phase while the profiles of the slower pulsars are 
not in phase ijKanba ch. 2002). Thus, a fast rotation rate seems to be instrumental 
in causing alignment of the pulses across the spectrum. 



1.2.3 Emission Models 

Particle acceleration inside the pulsar magnetosphere gives rise to pulsed 
non-thermal radiation and possibly also thermal emission from backflowing particles 
that heat the neutron-star surface. Rotating, magnetized neutron stars are natural 
unipolar inductors, generating huge vacuum electric fields. However, as was first 
noted by Goldreich fc JulianI l)l969|) in their classic paper, a rotating neutron star 
will not be surrounded by vacuum, since a large surface charge will build up on the 
star. This surface charge is unstable, because the induced vacuum electric field has 
a component parallel to the magnetic field at the stellar surface that exceeds the 
gravitational force by many orders of magnitude and is capable of pulling charges 
out of the star. If charges are not trapped in the surface by binding forces (i.e. if the 
crust is solid) a charge density, known as the Goldreich- Julian or corotation charge 
density, p^j ~ —uj ■ B/(27rc), builds up in the neutron-star magnetosphere. If the 
magnetospheric charge reaches the Goldreich- Julian value everywhere, it is able to 
short out the electric field parallel to the magnetic field (-En ), and the dipole magnetic 
field will corotate with the star. Corotation of the pulsar magnetosphere must break 
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Fig. 7.5. Multiwavelength pulse profiles llHarding et al.'. '2002b') of the Vela pulsar, 
as an example of emission characteristics of a middle-aged pulsar. 



down at large distances from the neutron star due to particle inertia. Where and how 
this happens is not yet understood, but it is certain that corotation cannot persist 
past the speed-of-light cylinder radius, Rlc = c/u), which is the distance at which 
the corotation velocity reaches c (see Figure 17^ . It is believed that in the outer 
parts of the magnetosphere there must be a transition to the wind zone, where the 
energy density of the particles is large enough to distort the poloidal dipole field into 
a toroidal relativistic wind flow that carries the spin-down energy of the pulsar, in the 
form of magnetic and particle energy, into a surrounding nebula wh ere (at least in 
the case of the Crab pulsar) it is dissipated as synchrotron radiation ( Rees fc GunnL 
1974HKennel fc CoronitiL ll984bV Although the standard picture of a pulsar magne- 
tosph ere described above is usually accepted on faith, global magnetospheric simula- 



tions ijKrause-Polstorff fc Micheilll98,4lp'etri et allEooilSpitkovskv fc AronsLl2002l) 

have not yet been able to show whether and how a pulsar magnetosphere reaches 
the nearly force-free (ideal MHD) state envisioned by Goldreich & Julian. However, 
if it is assumed that a for ce-free, corotating magnetosphere has been achieved, then 
Contopoulos et all l|l999(l have found a solution which smoothly connects the coro- 
tating magnetosphere to a relativistic wind. Inside the magnetosphere and at the 
light cylinder, the ratio (known as a) of the magnetic energy density to the particle 
energy density is large, primarily because according to the acceleration models to be 
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Fig. 7.6. Schematic 2D view of the high-energy emission geometry of several pro- 
posed models. 



discussed below, the particles receive only a fraction of the full potential drop across 
the open field lines. Therefore, a is large at the start of the wind flow, but must 
drop to 1 or below to accelerate the particles at the wind-termination shock. The 
pulsed emission is presumed to originate inside the corotating magnetosphere, and 
strong may develop to accelerate particles at two possible sites where E • B 7^ 0. 
These sites have given rise to two classes of high energy emission models: polar-cap 
models, where the acceleration and radiation occur near the magnetic poles in the 
inner magnetosphere, and outer-gap models, where these processes occur in the outer 
magnetosphere . 



7.2.3.1 Polar- C ap Model 

Polar-cap models toaugh ertv & HardingillQsilUsov fc MelroseLll99,4 l iDaughertv fc HardineL 

199fi^ advocate that particle acceleration occurs near the neutron star surface and 



that high-energy emission results from a curvature radiation or inverse Compton in- 
duced pair cascade in a strong magnetic field. There is some variation among polar- 
cap models, with the primary division being whether or not there is free emission 
of particles from the neutron-star surface. The subclass of polar-cap models based 
on free emission of particles of either sign, called space-charge limited flow models, 
assumes that the surface temperature of the neutron star (many of which have now 
been measured to be in the range T ^ 10^ — 10® K, see i l7.3(l exceeds the ion and elec- 
tron thermal emission temperatures. Although E^^ = at the neutron-star surface 
in these models, the space charge along open field lines abov e the surfa ce falls short 
of th e corotation charge, due to the curvature of the field jArons fc Scharlemannl 
19791) and to general relativistic inertial frame dragging ijMuslimo^^^Tsvgant Il99^ . 
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The E\\ generated by the charge deficit accelerates particles, which radiate inverse 
Compton (IC) photons (at particle Lorentz factors 7 ~ 10^ — 10®) and curvature 
(CR) photons (at Lorentz factors 7>10®). Because lower Lorentz factors are re- 
quired to produce pairs through IC emission, an IC pa ir formation front (PFF) wil l 
form first, close to the surface. However, it is found I Harding fc Muslimovl . 2002^ 



that the IC pair formation fronts do not produce sufficient pairs to screen the 
completely, thus allowing acceleration to the Lorentz factors ~ 10^ sufficient to pro- 
duce a CR PFF, where there are sufficient pairs to completely screen the £'||. The 
CR pair front will therefore limit the particle acceleration voltage and determine 
the high-energy emission luminosity. The CR photons can produce secondary elec- 
trons and positrons through one-photon pair production in the strong magnetic field 
near the neutron star surface. These pairs are produced in excited Landau states 
and radiate synchrotron photons which are energetic enough to produce more pairs. 
The combined curvature and synchrotron radiation from such pair cascades has been 
shown to produce high-energy pha se-resolved spectra sim i lar to observed spectra of 
X-ray and 7- ray Vela-like pulsars I Dauehertv Sz Hardind . 1996|) . The 7-ray spectra 



cut off very sharply (as a "super-exponential") due to one-photon pair production 
attenuation, at the pair escape energy, i.e. the highest energy at which photons 
emitted at a given location can escape the magnctosphcre without pair producing. 
A rough estimate of this cutoff energy, assuming emission along the polar-cap outer 
rim, ~ (2nr /cPy^"^ , at radius r, is 



2 GcVpi/2 (^)'^' maxjo.l, B,}, (^)'| (7.1) 

where i?o,i2 is the surface magnetic field in units of 10^^ G and R is the stellar 
radius. Although many high-energy photons are attenuated well below 1 GeV near 
the surface, in polar-cap cascades much of the radiation occurs at least 1-2 stellar radii 
above the surface, where the photon escape energy for pair production is several GeV 
(causing a spectral cutoff there) . At all but the highest fields there is a prediction that 
the spectral cutoff energy should be inversely proportional to surface field strength, or 
Bo = 6.4 X 10^9 {Ppyl'^ G for a dipole field at the pole. In fields above ~ 2 x lO^^ G, 
photon splitting, in which a single photon splits into two lower energy photons, 
becomes the dominant a ttenuation process and lowers the photon escape energy 
( Baring fc Hardind. I2OO1I) . The cascade radiation produces a hollow cone of emission 



around the magnetic pole, with opening angle determined by the polar cap half-angle, 
9 , at the radius of emission r. The polar cap opening angle is very small (a few 
degrees) unless the emission occurs more than a few stellar radii above the surface. 
But the pair cascades over most of the polar cap occur within several stellar radii 
of the stellar surface. Therefore, the wide double-peaked pulses of observed 7-ray 
pulsars like the Crab, Vela and Geminga, which require beam opening angles of the 
order of the magnetic inclination angle in this model, cannot be produced unless the 
pulsar is nearly aligned. 

Recently, this seemingly impossible requirement of polar-cap models has been 
eased somewhat by the realization that emission at the rim of the polar cap will 
occur at higher altitudes. The height of the CR pair front is within a stellar ra- 
dius of the neutron-star surface over most of the polar cap. However, near the edge 
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of the polar cap, at the last open field line which is assumed to be a conducting 
boundary, the parallel electric field drops to zero. Particles therefore need a longer 
distance to reach the Lorentz factor s neces s ary to produce pairs, causing the PFF to 
rise to higher altitude. A slot gap llAronat Il98i is formed as the PFF turns up to 
asymptotically approach the last open field line. The pair cascades that form on the 
edge of the slot gap at altitudes of several stellar radii above the neutron-star surface 
form a wide hollow cone of high-energy emission l|Mushmov fc Hardinellicioill . Even 
so, the inclination angles required in this model to reproduce the Crab and Vela 
profiles are less than 10° — 20°. This may be at odds with the ^ 60° viewin g angle 



inferred from the Chandra image of the Crab nebula I Weisskopf et al. . 2000(l . How- 



ever, acceleration of primary electrons in the slot gap can continue to high altitudes 
ijMuslimov fc HardineL l2004l) , where the emission may be viewed at larger angles due 
to the flaring of field lines. 

In polar-cap models, surface heating will occur as a result of the screening of the 
parallel electric field by pairs. Above the pair front, positrons will decelerate and 
some fraction will turn around and accelerate downward to the neutron-star sur- 
face. Because polar-cap heating emission is intimately tied to the pulsar acceleration 
mechanism, thermal X-ray emission provides very strong constraints on pulsar mod- 
els. In a space-charge limited flow model, the trapped positrons needed to screen the 
i^ll are only a small fraction of the number of primary electrons because this field 
was created by a small charge imbalance in the first place. The polar-cap heating 
luminosity is predicted by space-charge limited flow models to be only a few percent 
of the luminosit y of the non- thermal magnetospheric emission, and will fall in the 
soft X-ray band I AronsL 1981 1. Re cent self-consistent models of pola. r-cap heating in 
a space-charge limited flow model l|Harding fc Miislimo ^ 12001 [ 1200 have predicted 
heating fluxes from positrons trapped at CR and IC pair fronts. The heating from 
CR pair fronts predicts surface X-ray luminosity 



-{CR) 



10^^ ergs^i 



0.4 p-e/Vg-^/" iiP<Q.l B^/J, 



1.0 P-i/2 



if P>0.1 B, 



2' 
4/9 
0,12! 



(7.2) 



where tq is the pulsar characteristic age in units of 10^ yrs. Polar-cap heating 
may make the dominant contribution to thermal emission from old and middle-aged 
pulsars. The heating from IC pair fronts, with predicted X-ray luminosity 



-(IC) 



2.5 X 10^^ ergs"^ P"^/^ 



(7.3) 



can account for thermal emission components detected in the spectrum of some 
millisecond pulsars. 

Polar-cap models predict that emission from the open field line region is visible 
at all pulse phases. This is because the particles that are accelerated and initiate 
pair cascades at low altitudes radiate curvature (and possibly inverse Compton) 
emission at higher altitude on field lines extending to the light cylin der. Such off- 
pulse emission may have been detected a t 7-ray energie s by E GRET <|Fierro et all 
1998(l and at X-ray energies by Chandra ( Tennant et all 2001 ) . 
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7. 2. 3. 2 Outer- Ga y Model 

Outer-gap models 1 Cheng et al. . 1986t RomaniL 1996^ assume that acceler- 
ation occurs in vacuum gaps that develop in the outer magnetosphere, along the 
last open field line above the null charge surfaces, where the Goldreich-Julian charge 
density changes sign (see Figure 17.6(1 , and that high-energy emission results from 
photon-photon pair production-induced cascades. The gaps arise because charges 
escaping through the light cylinder along open field lines above the null charge sur- 
face cannot be replenished from below. Pairs from the polar-cap cascades, which 
flow out along the open field lines, may pollute the outer gaps to some extent (and 
vice versa), but this effect has yet to be investigated. The electron-positron pairs 
needed to provide the current in the outer gaps are produced by photon-photon pair 
production. In young Crab-like pulsars, the pairs are produced by CR photons from 
the primary particles interacting with non-thermal synchrotron X-rays from the same 
pairs. In older Vela-like pulsars, where non-thermal X-ray emission is much lower, 
the pairs are assumed to come from the interaction of primary particles with thermal 
X-rays from the neutron-star surface. Some of the accelerated pairs flow downward 
to heat the surface and main tain the required thermal X-ray emissio n. The modern 
outer-gap Vela-type models ( Romani 1996t Zhang k, CheneL 1997|) all adopt this 
picture. 

Although there seems to be agreement on the radiation processes involved in the 
outer gap, the full geometry of the gap is still not solved. Two approaches to such a 
solution are currently underway, but nei ther is near to defining the complete three- 
dimensional gap geometry. One group I Zhang fc Chenel 1997t Cheng et al. . 20001) 
solves the ID Poisson equation perpendicular to the magnetic field lines, result- 
ing in a gap geometry for young pulsars that is a long, thin sheet bounded be- 
low by the last open field l ine and above by pro duction of pairs. The other group 
I Hirotani fc ShibataL 2001 : Hirotani et all 2003(l obtains solutions to the ID Pois- 
son equation along the magnetic field (assuming a gap width across field lines) and 
finds that the gap is limited parallel to the field by pair creation. The actual gap 
ge ometry is probably sornewhere in between. The long n arrow outer-gap geometry 
of lRomani & Yadigar and lCheng et alJ has been successful in re- 

producing the observed double-peaked pulse profiles for large inclination angles. In 
these models, the electrons/positrons flowing outward/inward along the last open 
field line radiate outward/inward from the null surface. The emission forms a wide 
fan beam which produces in general a wide double-peaked profile, much like the 
observed Crab profile. The narrow peaks are the result of caustics formed in the 
observer's frame by the cancelation of the phase shifts from relativistic aberration 
and light travel time, so that emission from a large range of altitudes is compressed 
into a narrow phase range. However, due to the fact that no outward emission origi- 
nates below the null surface, the profile falls off very abruptly at the outer edges (see 
Figure W7l\ with no leading and trailing or off-pulse emission. The inability to pro- 
duce trailing and off-pulse emission, as is seen in high-energy p ulsar profiles, is thus 
a serio us p roblem for curr e nt ou ter-gap models. The models of Hirotani fc ShibatU 



1 200l[l and iHirotani et al.l l|2003f) allow external currents to fiow through the outer 
gaps, thus producing a possible extension of the gap below the null surface. In the 
case of the Crab however, the extent and magnitude of the emission below the null 
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surface is small and would not produce much ofF-pulse emission. The formation of 
the leading peak in double-peaked profiles is also problematic since it requires emis- 
sion very close to the light cylinder, where the structure of the magnetic field has 
not been determined. 

The spectrum of radiation from the outer gap is a combination of synchro-curvature 
radiation from pairs inside and outside of the gap and inverse Compton scattering of 
the synchrotron emission by the pairs. Unlike in polar-cap models, pair production 
in outer-gap models is essential to the production of the high-energy emission: it 
allows the current to flow and particle acceleration to take place in the gap. Beyond a 
death line in period-magnetic field space, and well before the traditional radio-pulsar 
death line, pairs cannot close the outer gap and the pulsar can not emit high-enerCT 



radia tion. This outer-gap death line for high-energy pulsars ijChen fc RAidermanL 
falls around P = 0.3 s for B ^ 10^^ G. Geminga is very close t o the outer-gap 
death line and recent self-consistent models l|Hirotani et al. have difficulty 



accounting for GeV 7 rays from pulsars of this age. Polar-cap models, on the other 
hand, predict that all pulsars are capable of high-energy emission at some level, so 
that detection as a radio-loud high-energy pulsar is thus a matter of sensitivity. 

Half of all the particles produced in the outer gaps are accelerated back toward 
the neutron-star surface and radiate curvature photons that can produce pairs in 
the s t rong magnet i c field near the surface, initiating pair cascades I Zhang fc Chend . 



1997t Wang et al. . 1^98*). But the pairs still have enough residual energy to heat 



the surface at the footpoints of open field lines that thread the outer gaps. The 
resulting thermal emission has T ^ 1 keV, which is much higher than the thermal 
emission temperatures observed in pulsars. According to the model, this emission 
is not observed directly (except right along the poles) but only through the blanket 
of pairs produced by the downward-going particle cascades. The 1-keV photons 
are reflect ed back to the surface by th e pair blanket through cyclotron resonance 
scattering i Haloern fc RudermanL 1993j) . and are re-radiated from the entire surface 



at a temperature around 0.1 keV. Thus, these outer-gap models predict three X-ray 
emission components: hard thermal emission from direct heating of the polar caps 
(seen only along the poles), soft thermal emission reflected from the pair blanket, 
and non-thermal emission from the downward pair cascades (Cheng & Zhang, 1999). 
The components actually observed from a particular pulsar depend on inclination 
and viewing angle. The predicted X-ray luminosities for pulsars in the ROS AT hand 
can account for the observed Lx = 10~^E relation. 



7.2.4 It's The Geometry... 

Theoreticians have spent years building pulsar emission models by start- 
ing from the fundamental electrodynamics of the particle acceleration and radiation 
processes, predicting observable pulsar characteristics from the bottom up. This has 
had some success, but both polar-cap and outer-gap models as currently formulated 
also have some basic, unresolved problems. Ultimately, it is the emission geometry 
required by the observations, which tells us the distribution of radiating particles, 
that will drive us to the correct understanding of pulsar high-energy emission. 

Putting together the results of a number of attempts, in different types of models, 
to reproduce the narrow and often double-peaked pulse profiles of high-energy pul- 
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sars, one common thread seems to emerge. Nearly all successful geometrical models 
of pulsar high-energy emission have assumed enhanced particle acceleration and radi- 
ation a long the last open field lines o f a magnetic dipole. In their polar-cap emission 
model, Daughertv fc Hardind ljl99(](l were led to assume strongly enhanced particle 



flow on the rim of the polar cap. Acceleration in the slot gap along the polar-cap 
rim later provided a physical basis for this i dea |Mus l imov fc Ifard ing, 2003). The 
outer-gap models of Romani fc Yadigaroelul ( 1995(l and lCheng et al.l |2000|l assumed 



acceleration and emission essentially only along the last open field line between the 
null charge surface and the light cylinder in a thin outer gap. But these models can- 
not produce the non-thermal off-pulse emission seen through the who le pulse phase of 
the Cr ab optical and Vela and Geminga 7-ray profiles. Very recently, Dvks fc Rudald 



l|200,'y have explored a purely geometrical emission model, known as the "two-pole 
caustic" model in which they assume that particles radiate all along the last open 
field lines of both poles, from the neutron-star surface to the light cylinder (see Fig- 
ure l7.6|l . Taking into account relativistic effects of aberration and light travel time, 
the phase shifts of photons emitted at different altitudes along the trailing edge of 
the polar open field line region cancel, so that the emission in the observers frame 
is compressed into a small range of phase, forming a caustic peak in the profile. 
The caustic peaks in the two-pole caustic model have the same origin as the caustic 
peaks that form from outward emission above the null charge surface in the outer- 
gap models. What makes the two-pole caustic model fundamentally different is the 
addition of outward emission below the null surface, that allows an observer to view 
caustic emission from both poles as the neutron star rotates past, in contrast to 
outer-gap models where an observer can view high-energy emission from only one 
pole. The resulting profiles for viewing angles crossing caustics from both poles are 
astonishingly similar to the observed profile of the Crab (see Figure I7T|) and several 
other high-energy pulsars, and seems to do better than most of the other 'bottom-up' 
models in reproducing the observed geometry. This seems to be telling us that both 
conventional polar-cap and outer-gap models are missing something crit i cal in their 
description of the basic electrodynamics. Recentlv. lMiislimov fc Hardin have 



shown that high-altitude radiation from the slot gap could form a physical basis for 
the two-pole caustic model, producing very similar profiles. 

The geometry of emission at different wavelengths is also different in various mod- 
els. In polar-cap models, the high-energy and radio emission are physically con- 
nected, since the electron-positron pairs from the polar-cap cascades are thought to 
be a necessary ingredient for coherent radio radiation. However, the relative geom- 
etry of the two emission re gions is not very constrained bv observati on. Models for 
radio emission morphology ( RankinL 199,'^ Ilvuc fc Manchester! Il988l) consist of core 



and conal components, emitted within a distance of about 10-100 stellar radii from 
the surface. If the polar-cap high-energy emission cone occurs at a lower altitude 
than the radio emission, as would be expected, then the leading edge of the radio cone 
would lead the first high-energy peak. The trailing edge of the radio conal emission 
would then have to be undetected in the CGRO p ulsars. This type of radio geometry 



for young pulsars was suggested bv iManchesten (1996) and a class of pulsa rs with 



one-sided radio conal emission was proposed by iLvne fc Manchester! l|l988l) . This 
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pict ure has recently received some observational support by ICrawford et al. I il2nnih 
and lCrawford fc KeimI ^00^. 



High-energy emission in the outer gap is generally radiated in a different direction 
from the radio emission, which allows these models to account for the observed 
phase offsets of the radio and high-energy pulses. At the same time, there will be 
fewer radio-high-energ y coincidenc es and thus a larger number of radio-quiet high- 



energy pulsars. In the Romani fc Yadigaroglul |l995() geometrical outer-gap model. 



the observed radio emission originates from the magnetic pole opposite to the one 
connected to the visible outer gap. Many observer lines-of-sight miss the radio beam 
but intersect the outer-gap high-energy beam, having a much larger emission solid 
angle. When the line-of-sight does intersect both, the radio pulse leads the high- 
energy pulse, as is observed in mos t T-ray and soni e X-ra y pulsars. 

In the two-pole caustic model of Dvks fc RiidalJ l|200,ll) . trailing-edge radio conal 



emission would arrive in phase with the high-energy emission due to the fact that 
emission from a wide range of altitudes is compressed into a narrow phase range to 
form the caustic. This may be a nice explanation for why the radio and high-energy 
emission tends to be phase-aligned in some of the fastest rotators, such as the Crab 
and the millisecond pulsars. 

7.2.5 Polarization Properties 

Phase-resolved polarimetry of rotation-powered pulsars has had enormous 
diagnostic capability at radio and optical wavelengths and could also be a powerful 
diagnostic in the X-ray range. Several X-ray polarimeters are planned for the near 
future (AXP, POGO and MEGA) which could measure X-ray polarization charac- 
teristics of the brightest sources, such as the Crab pulsar, Cyg X-1 and Her X-1. 
The pulsed non-thermal radiation from relativistic particles in the magnetosphere 
is tightly beamed along the neutron-star magnetic field lines and thus the emitted 
radiation is believed to be highly polarized either parallel or perpendicular to the 
field lines. Since the field well inside the speed of light cylinder rotates as a solid 
body with the star, measurement of the polarization properties as a function of pulse 
phase can provide a multidimensional mapping of the pulsar emission. The expected 
signature of emission near the poles of a dipole fiel d, an 'S'-shaped swing of the po - 
larization position angle through the pulse profile I Radhakrishnan & CookeL 19691) . 



has been seen from many radio pulsars and has generally been taken as proof that 
the radio emission originates from the open field lines of a magnetic dipole. 

The different high-energy emission models share the common emission mechanisms 
of curvature, synchrotron and inverse Compton radiation from highly relativistic par- 
ticles, albeit in very different locations. These mechanisms all intrinsically produce 
highly polarized radiation (up to 70%, depending on the particle spectrum). Each of 
the models predicts a distinct variation of position angle and degree of polarization 
across the pulse profile. Figure FTTI shows simulations of expected position angle and 
percent polarization as a function of pulse phase for outer-gap and polar-cap mod- 
els of the Crab pulsar . Simulations are also shown for the two-pole caustic model 
I Dvks fc R,uda5l200.'^ . The predicted position-angle swing in the polar cap model 



shows the classic S-shaped pattern of the rotating vector model since the observer is 
sweeping near the magnetic pole (phase 0). The most rapid change in position angle 
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Fig. 7.7. Pulse profiles, polarization position angle and degree of polarization as 
a function of pulse phase for the Crab pulsar. The left panel the observe d optical 
pulse profile and polarization characteristics as measured by OPTIMA JKellneA 
l2002f) . The right-hand panels are predicted profiles and corresponding polarization 
characteristics of t he polar-cap, two-pole caustic and outer-gap model geometries 
JDvks et alJ.l200il . 



occurs near the magnetic axis and there is no predicted decrease in degree of polar- 
ization. The two-pole caustic model predicts rapid swings of position angle through 
the peaks, due to the large change in field orientation, as the viewing angle sweeps 
past the pole (phase 0) and approache s the trailing edge field line that is the origin 
of the caustic peaks l|Dvks et al.Ll2004|) . There is also a predicted decrease in percent 
polarization on the trailing edges of the peaks, as the observer is viewing emission 
in the caustics and from rotationally distorte d field lines near the ligh t cylinder that 
fold back and overlap. The outer-gap model i Romani fc YadigarogluLil995|) predicts 
position angle swings through the peaks, though they are smaller and occur only for 
carefully selected combinations of inclination and viewing angle ( Dyks„e t al., 2004). 
There would also be a predicted dip in degree of polarization, now primarily at the 
second peak, since the second peak in the outer-gap model originates from the same 
caustic that is the origin of the dips at both peaks in the two-pole caustic model. 
The signature of caustic emission is thus a dip in the percent polarization and a rapid 
swing of the p osition angle at the pulse peak. A recent measurement by OPTIMA 
of optical polarization position angle and degree of polarization as a 
function of pulse phase for the Crab pulsar exhibits sharp swings in position angle in 
the peaks and dips in degree in polarization just following each peak. There appears 
to be a constant, unpulsed background component (at about 120° position angle and 
30%) which if subtracted from the data, best matches the predictions of the two-pole 
caustic model. 
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7.3 Thermal Emission from Cooling Neutron Stars 



7.3.1 Motivation and Background 

The existence of magnetospheric emission from neutron stars, as described 
in i )7.2l came as an observational surprise to neutron-sta r theorists. By contrast, not 
long after the existence of neutron stars was proposed I Baade fc Zwickvl 19341) . it 
was also pointed out that such a star would form hot, with surface te mperatures over 
10^ K, and that surface emission might be detectable I Zwickvl 19381) . The discovery 
of cosmic X-ray sources in the 1960's provided hope that neutron stars might one 
day be observable, and detailed mode li ng of their po tential observational properties 



began l|Mortonl [l96i iTsurutal 11964 - iBahcall fc W oU (1965J pointed out that the 



equation of state of matter at densities comparable to those in atomic nuclei, which 
cannot be studied in terrestrial laboratories, might be constrained by observations 
of neutron-star cooling. This is because the magnitude of the stellar neutrino lumi- 
nosity, which, in the early history of the neutron star, should far exceed that of the 
photon luminosity, should depend strongly on what species are present in the core, 
and on how they interact. Prodigious neutrino emission would result in more rapid 
cooling than photon cooling alone. Studying neutron-star cooli ng began in quantati- 
tive earnest with the computation of the first cooling curves ( Tsuruta fc CameronL 
196^. 



The first actual detections of thermal emission from the surfaces of neutron stars 



came much later, fo l lowing the launch of the Ei nsteAn llHelfancj ^\ ^l.|,|l98Cl|:|Ha rnden et al 
19851 iMatsiii et all . Il988l) and EXOSAT (e.g. Brin kmann fc Qgelmanl 119811) X-ra y 
observatories in the 1980's. The situation pre-1990's was reviewed bv lTsurut arill98(il) 



The launch of ROSATva 1990 heralded the next major progress (e.g. iFinlev et al 



1992H6gelman et al.l [l99i iHalpern fc Ruderm'^ Il993l). with mport ant contribu- 
tions made by ASCA (e.g. iHalpern fc Wanel Il997l IZavhn et all Il998|) as weU. By 
1997, four bona fide detections of thermal emission from young or middle-aged pulsars 
had been reported (from PSRs B0656-K14, B08 33-45, B1055-52, and th e Geminga 
pulsar) along with interesting upper limits (see Becker Triimpeil. Il99^ and refer- 
ences therein) . The recent launches of the C handra and XMM-Newton observatories 
have already yielded important results (e.g. IPavloy et all boOlt iMarshall fc Schul j . 

rinnollT^ ij. i _l I lonnrJIo 1 _i _ l ' 'r>r>nr>i ' 'n: : .j- -l''nnnnl n j. 



2002HKargaltsev et alll2002HSanw al et aL 2002b': 'Bignami et al.','20 0d). Recent re - 
views of relevant observ ations include . Becker fc Pavlov, ( 2002.') : .Pavlov et alJ ll2002bl) : 
Pavlov fc Zavih] tooi). Recent reviews of theory include iTsurutal l|l998|) . IPagd 



19981). andlYakovlev et al 



1 19991) . 



As discussed in ^7.2\ thermal emission can also be produced as a byproduct of 
emission processes in the magnetosphere. Such emission is believed to be unrelated 
to that originating from initial cooling. Other heating mechanisms are predicted 



to rnatter only in the photon-cooling era, i.e. for t > 10 yr (e.g . lAlpar et al 
19841 lUmeda et all 11993: Ivan Riper et allll995l:lLarson fc LinkLll999|) 



One notable 

exception are magnetars, in which the decay of the field is thought to provide a major 
energy source that dominates over all others (see chapter 14). 
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7.3.2 Theory of Neutron Star Cooling 

Neutron stars are born in the violent collapse of stars having ZAMS mass 
~8-20 Mq, which have core temperatures prior to collapse of Tc ^ 10^ K. Following 
collapse, it is expected that the neutron star forms with central temperature Tc ~ 
10"'^"'^ K, although prodigious neutrino emission results in Tc ~ 10^ — lO^*' K within a 
day or so. Neutrino emission dominates over photon emission until Tc ^ 10^ — 10® K 
(depending on leading neutrino emission process), which is reached at least ~ 10'^ yr, 
and more likely ~ 10^ — 10^ yr after birth. The thermal behavior before 10^ yr 
depends strongly on cooling model. After t ~ 10^ yr, various internal and external 
heating mechanisms dominate thermal emission over the initial cooling. Therefore, 
young, i.e. having age S; 10^ yr, neutron stars for which the surface temperatures 
are expected to drop either slowly or quickly from 2 x 10^ to ^ 2 x 10^ K, with 
core temperatures dropping from Tc ^ 10^ to 10^ — lO'^ K, are interesting targets 
for studies of neutron-star interior physics. Note that in spite of these apparently 
high temperatures, from a statistical mechanics point of view, neutron stars can be 
considered to be at T = 0, since for their average densities, the Fermi energy is 
~ 100 — 500 MeV, while the internal temperature of a one-year old star is certainly 
below 0.1 MeV. Next we consider in turn the main elements that determine how a 
neutron star cools. 



7.3.2.1 The Stellar Core 

During the first ~ 10^ — 10^ yr, a neutron star cools mainly via neutrino 
emission from its innermost regions (which have densities p ^ 10^" g cm^"^). For ages 
t 10 — 50 yrs, the internal layers are isothermal. They have an enormous density 
gradient which results in different neutrino emissivities at different radii. As a rule, 
the most efficient neutrino emission is produced in a stellar core, which extends from 
p « po/2 to the stellar center p ^ (3 — 15) po. Here, po « 2.8 x 10^"* g cm~'^ is the 
density of nuclear matter in atomic nuclei. The composition and equation of state 
of supranuclear matter (p > 2 po) in neutron-star cores are largely unknown but are 
of considerable interest to basic physics. 

The main cooling mechanism in the core is neutrino emission via the "Urea" pro- 
cess*. In its simplest version adopted for neutron-star cores, the Urea process is a 
sequence of a /?-decay and an electron ca pture, n — > p + e+V^ , and p + e — > n + i>c- 



This is the so-called direct Urea process l|Lattimer et al.Ul991l) . However, to proceed 



at interesting rates in the neutron-star interior for a wide range of expected interior 
properties, the Urea process requires a spectator reactant to simultaneously conserve 
energy and m omentum. The proces s with the spectator is called the "modified" 



p roces s 

Urea process (jChiu fc Salpeten 19641) . Neutrons, protons, and other baryon spec- 



tators greatly suppress the direct Urea rate (typically, by 5-7 orders of magnitude) 
because of the lack of phase space due to degeneracy. The most efficient direct Urea 
cooling, without any spectator, is possible if the proton fraction is sufficiently high. 
This may happen in the cores of sufhcicntly massive stars. If hyperons are pre sent 
in the core, they initiate additional direct Urea processes i Prakash et al. . 1992|) . as 



Urea is the name of a long-since-closed casino in Rio de Janeiro, and was adopted as a name for 
these reactions by Gamow & Schoenberg (1941) who saw a parallel between how casinos extract 
money from players and how nature extracts energy in these reactions. 
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efficient (or nearly as efficient) as the basic direct Urea process with nucleons. If the 
direct Urea process is forbidden (e.g., if the proton fraction is insufficiently high) but 
the dense matter contains pion or kaon condensates, weaker processes similar to the 
direct Urea on es (though stronger than the modified Urea ones) can open (see, e.g. 
Pethickl I992I for a review). 



7. 3. 2. 2 Equation of State 

An equation of state (EOS) is the pressure-density-temperature relation- 
ship of matter. In a neutron star, most of the mass is at densities two or three 
times pq (the nuclear matter density). At these densities, the EOS is unknown. 
Nuclear physics can constrain some aspects of the high-density EOS, however cer- 
tain fundamental parameters, such as the compression modulus, th e bulk syrnmetr'y 



tam lunaamental parameters, sucn as tne compression modulus, th e bulk syrnmetry 
energy and the effective n ucleon mass remain weakly constrained I Lattimeil . 19921 : 



Lattimer fc PrakashI 2001 ) The EOS is crucial to understanding neutron-star struc- 



ture; it determines, among other things, the mass-to-radius relationship, the stellar 
binding energy, the stellar moment of inertia and the relative moments of different 
components, the minimum and maximum masses, and the maximum angular veloc- 
ity. From a cooling point of view, the most important parameters determined by 
the EOS are the density, pressure and temperature profiles of the star, as well as its 
composition. The profiles are important because, for example, the location, amount 
and type of superfiuidities present depend on the density profile, via a critical tem- 
perature/density relation (see i l7.3.2.3|l . The composition is important as different 
species of particles will react differently. 

Many EOSs have been proposed, and the range of possible properties of neutron- 
star matter is large. For example, the range of predicted pressure of matter having 
density twice po is approximately a factor of five, depending on the choice of plausible 
EOS, where plausible means consistent with available laboratory data on atomic 
nuclei and nucleon scattering ( Lattimer fc PrakashI 200lj) . Phase transitions to kaon. 



pion or hyperon-rich matter are all possible depending on the core d ensity. 



Fig ure rr^ fleftl shows predicted cooling c urves for a commonly fe.g. lYakovlev et al 



l200l[l assumed core EOS (originally from iPrakash et al.L 11988^ for three different 
masses. Curves for M < 1.35 Mq are very similar to that for M = 1.35 Mq. Above 
this mass (to be more exact, above 1.358 Af©), for this EOS, direct Urea cooling sets 
in; this is clear in the curves. Also shown in the Figure are two schematic curves 
assuming a pion condensate (but no direct Urea if the pion condensate were absent) 
in the core at Af > 1.35 Mq. The effect is clearly enhanced cooling, but not as much 
as for direct Urea. 

7.3.2.3 Superfluidity and Superconductivity 

Proton and neutron superfluidity in the neutron-star interior can play a very 
significant role in cooling. (Note that proton superfluidity means superconductivity.) 
This is because it suppresses traditional neutrino processes in volving nucleons. On 
the other hand, it produces a specific type of neutrino emission | Flowers et al. , 1976(l 



associated with Cooper pairing of nucleons. Superfluidity also affects nucleon heat 
capacity. However, its onset, determined by a critical temperature, is very poorly 
known at neutron-star densities. Overall, superfluidity in different forms can mimic 
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Fig. 7.8. Neutron-star cooling curves for a variet y of parameters and rn odels (cour- 
tesy D. Yakovlev). All curv es are for EOS "A" of lYakovlev et alJ i2001[) . or "model 
I" of IPrakash etai] itTosl) with compression modulus of saturated nuclear mat- 
ter 240 MeV. Left: Curves for non-superfluid stars for three masses (labelled in 
solar masses) are shown with solid lines. Dotted lines are for pion condensate in 
the core, for 1.4 and 1.7 Mq stars (see H7.3.2.2II . Middle: Curves for superfluid 
neutron stars of masses 1.35 and 1.5 Mq. Four curves are shown for each mass: 
solid curves are for non-superfluid sta rs (as in the lef t plot) ; dotted curves are 
for proton superfluidity model "Ip" of iKaminker et all i2002f) : dashed curves are 
for proton superfluidity model "2p" of the same reference; dot-dashed curves are 
for the combination of proton superfluidity model "pi" and neutron superfluidity 
model "3nt," also from Kaminkcr ct al. (20Q3)- Note that the Ip and 2p curves for 
the 1.35 Mq star coincide as both suppress modified Urea cooling. Also, note that 
the 2p curve for the 1.5 AIq star sits among the 1.35 Mq star curves. Right: Curves 
for non-superfluid stars for masses 1.35 and 1.4 Mq for different accreted envelope 
masses and dipolar magnetic fields: solid curves are for zero accreted mass and zero 
magnetic field B (as in the left plot); dotted curves are for B = 10^* G at the pole; 
dashed curves are for B = and accreted envelope mass 10"'^ Mq. The dot-dashed 
curve is for a 1.35 Mq solar mass star having 5 = and accreted envelope mass 
10"'^ Mq, but with proton superfluidity included. It is shown to demonstrate the 
approximate upper bound for temperatures of neutron stars having ages between 
~1 and ~100 kyr. 



enh anced cooling, challeng ing the promise of using cooling to determine the EOS 



See Kaminker et al. ( 2002() and Yakovlev et al. 12004^ for recent reviews of the 



subject. 

The relevance of superfluidity can be summarized as follows. 



In very low-mass 
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stars (M < 1.35 Mq for the EOS assumed here), the effect of core singlet-state ("^So) 
proton superfluidity is well determined and nearly independent of mass or EOS. 
These stars cool very slowly as even modified Urea processes are suppressed, so that 
nucleon-nucleon bremsstrahlung, a slower process, dominates. The hottest young 
neutron stars are therefore plausibly explained as having the lowest mass. On the 
other hand, for very high-mass neutron stars {M ^ 1.5 Mq), proton superfluidity 
becomes inefficient, and one expects very rapid, model-independent cooling, because 
direct Urea is too strong to be suppressed. A neutron star that has cooled so rapidly 
has not yet been observed (see H7.3.4(l . 

For intermediate-mass neutron stars, the situation is highly model dependent, and 
by assuming an EOS and proton superfluidity model, one can, at least in prin- 
ciple, use observed surface temperatures to determine neutron-star masses (e.g. 
Yakovlev et all 2fl0llh . In these stars, direct Urea processes would proceed unim- 



peded were it not for moderate proton superfluidity supp ression. Figur e lT.SI fmiddle) 



shows two different models for proton superfluidity (see iKaminker et al.^ .2002. and 
references therein for definitions). Mild pairing of neutrons in the stellar core 
results in very strongly enhanced cooling of stars with M < 1.35 fFigure I7!H1 
middle) that is inconsistent with the observations thus far (see i )7.3.3|l . 

7.3.2.4 The Stellar Envelope: Composition and Magnetic Field 
Interpreting observed thermal emission from neutron stars means mapping 
the observed surface temperature Ts to that of the stellar core Tc. The temperature 
gradient from the surface inward is mainly in the stellar envelope, the layer beneath 
the atmosphere down to the isothermal internal region which effectively acts as a 
thermal insulator for the bulk of the star. The insulating envelope is commonly 
defined as the region having density p < 10^° g cm^'^. This has thickness only a 
few tens of metres, much smaller than the stellar radius. The insulating envelope 
thickness decreases as the star cools. 

For envelopes made of iron, for T, ~ 10^ K, T , ~ 10^ K with T,. 2 T^^ for a 



typical surface gravity l Oudmundsson et all Il983lj . How ever, the mapping is highly 
sensitive to envelope composition. IChabrier et al. show that even a tiny (~ 



10""'^'^ Mq) amount of accreted light-element matter, such as from fallback after the 
supernova, from the interstellar medium or from a binary companion, can have a 
substantial effect on Ts for a given Tc. Specifically, accreted envelopes have lower 
thermal conductivity. As a consequence, they lead to higher Ts in the neutrino- 
cooling era but lower Tg in the photon-cooling era. Figure lTTHl fright) shows the effect 
of an accreted envelope on predicted cooling curves of two different mass neutron 
stars. 

The magnetic field in the envelope can also plav a major, though complicated, role 
' Ri ner fc T.ambl.ll 98lHPag4 |l flfl5':'Shib anov Y akovW. h flflfilffT^v] Hernnuistl 



11997; ■Potekhin fc Yakovlevli20nillHevl fc Hernauistl . .200Ll . A magnetic field affects 
all plasma components, especially the electrons. In particular, electron motion per- 
pendicular to the field lines becomes suppressed by classical Larmor rotation. In 
addition, this motion is quantized in Landau levels. Classical and quantum effects 
greatly modify the electron thermal and electrical conductivities. Thus the effect of a 
magnetic field on the thermal structure of the insulating envelope depends strongly 
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on the field geometry, with radial fields effectively reducing the insulation by the 
envelope (quantum effect), and fields tangential to the surface increasing it (clas- 
sical effect). Thus a dip olar field will have a very d ifferent effect than an (often 
assumed) radial field fsee IPotekhin &: Yakovlevl 2001 , for a review), with the tan- 
gential and radial regions potentially cancelling local deviations when averaging over 
the surface (P age, 1995; Shibanov & Yakovlev, 1996; Potckhin & Yakovlcv, 200ll 
[Hevl fc Hernouist). boOli: jPotekhin et all |200^. Figure O (right) shows the effect 
of a dipolar magnetic field of 10"'^'* G (at the pole) on cooling curves for two different 
masses. Curves for 10^^ G are close to the zero field case. The net effect is clearly 
much more subtle than for light elements, mass, or EOS. Note that the above consid- 
eration assumes to active magnetic field decay of the type invoked in the "magnetar" 
model (see Chapter 14). 



7.5.^.5 The Stellar Atmosphere 

The neutron-star atmosphere, defined as the region having density p ;5 
10^ g cm^'^, is typically only 0.1-10 cm thick, yet ultimately determines many of 
the properties of the emerging thermal photon flux. Although in studying thermal 
radiation the assumption of a blackbody spectrum is often made, every realistic at- 
mosphere model predicts a significant deviation from a Planckian curve. The nature 
and degree of the deviation, both from the continuum and in the form of lines (see 
i i7.3.3.4|l depends strongly on atmosphe ric composition a. s well as on magnetic field. 
A nice review of this subject is given by Zavlin fc Pavlov! 1 20021) . We summarize the 
basic issues below. 

An unadulterated neutron star might be expected to have an iron atmosphere, 
in keeping with the pre-coUapse composition of the progenitor core. Heavy-element 
spectra, with even tiny (i.e. ~0.05%) metallicities, are expected to contain many 
spectra l features and photoion ization edges produced by ions at various stages of ion- 
ization ( Romani 1987[.RaiagODal fc Romani 1996t Zavlin et al. , 1996t Raiaeopal et al 



1993). However, blurred to lower resolution, a heavy-element atmosphere resembles 
a blackbody of comparable effective temperature. Figure 17.91 shows sample atmo- 
spheres containing different metal abundances, for comparison with a blackbody 
spectrum. 

On the other hand, a spectrum emerging from a light-element (H or He) atmo- 



spher e is predicted to have few absorption lines at X-ray energies (e.g. iHo et al 
2003|) - in fact none if the atoms are fully ionized, apart from cyclotron lines. Nev- 



ertheless, such an atmosphere has a clear signature: a high-energy tail in the spec- 
trum, resulting in a significant overestimate, by as much as a factor of ^ 3, of the 
surface temperature if a blackbody model is fitted instead (see Fig. 17.91 RomaniL 
1987t Raiagopal fc RomaniL ll99(iHZavlin et all Il99(il) . For example, fitting a black- 
body model to the X-ray spectrum of the V ela pulsar yields kT^ = 127 eV, versus 

. Also, blackbody 



59 eV for a hydrogen atmosphere l)Pavlov et all |2 



fits to hydrogen atmosphere spectra yield emitting radii that are much smaller from 
those inferred with hydrogen atmosphere models; for example, Marshall fc Schui3 
l|2002() . as revised bv lBrisken et all l|2003b(l . find = 156 km and 8.5 km for hy- 
drogen and blackbody fits to the X-ray spectrum of PSR B0656-I-14, respectively. 
The high-energy tail is produced because light-element opacities decrease rapidly 
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Fig. 7.9. Spectra of emergent radiation in non-magnetic neutron star atmospheres 
having Tcff = 10® K for different metalUcities Z {Z — 2.0% corresponds to solar 
metallicity.) The corresponding blackbody is shown with a dotted hne. From 
fZavlin fc Pavlov (2002). 

with increasing energy, so that at higher energies one is observing deeper layers of 
the star which are at higher temperatures. Although an unadulterated neutron star 
should have an iron atmosphere, even tiny amounts (e.g. 10~^° -^0) of accreted 
interstellar material should result in a pure hydrogen atmosphere, because of the 
short time scale for gravitational settling. 

As in the envelope, an appreciable magnetic field can also have a major influence 
on the atmosphere. The most important issue is the huge difference that exists be- 
tween the opacities for different polarizations. Specifically, the magnetic atmosphere 
is much more transparent to photons with electric-field vector perpendicular to the 
k ■ B plane than is a non-magnetic atmosphere. Thus, the emergent radiation is ex- 
pected to be highly polarized (see tl7.^{.3.3(l . Furthermore, the emergent flux depends 
on the direction of the magnetic field; hence for non-radial field geometries pulsations 
are expected even for uniform temperature (which is, in any case, implausible, since 
variations in surface magnetic field strength result in temperature variations). 

For magnetic field strengths greater than a quantum critical field Bq = 27rTOgC'^/ (eh) = 
4.414 X 10"'^^ G, so-called "vacuum polarization," that is the effective dielectric con- 
stant of the vacuum, must be accounted for in the radiative transfer equations, given 
the plasma densities in the envelope. This is interesting physically as, in spite of 
being predicted by QED theory, this effect has not been confirmed experimentally. 
The net effect of vacuum polarization is to soften the high-energy tail of the hydrogen 
spectrum, though the emission is still harder than the blackbody spectrum of the 
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same effective temperature | Ho fc LalboOlilOzeA 2001 ; Zane et al. , 2001 ; Lai fc HoL 



20021 IHo fc Lai l20n.^: Ic^^ell l2003l iLai fc Hoi. 1200.4 .' 



7.3.3 Determinations of Neutron- Star Temperatures and Ages 

7.3.3.1 Measuring Neutron-Star Temperatures 

There are several challenges to measuring neutron-star temperatures. First, 
the temperature measured by a distant observer T^ is not that at the neutron 
star surface Ts. This is due to gravitational redshifting; the two are related via 
J.OO _ Ts[l — 2GM/{RcP)Y/^. Here, i? is a circumferential neutron-star radius (which 
determines the proper length 2t:R of the stellar equator) . Similarly, the stellar lumi- 
nosity as measured by a distant observer is given by — 4:TTaT^R'^[l — 2GM/{Rc'^)], 
where the 7 specifies that this is the bolometric photon (as opposed to neutrino) lu- 
minosity. The apparent radius which would be measured by a distant observer if the 
latter could resolve the star is R°° = R/[l-2GM/{Rc^)]^^^. Thus the determination 
of Ts from T^ requires knowledge of M / R, which is unknown, and which depends 
on the unknown EOS in the stellar core for any M or R that could be determined 
independently. 

Particularly problematic in the pre- Chandra era is that young neutron stars that 
are pulsars with high E are ubiquitously associated with pulsar wind nebulae (see 
ii7.6|l . These extended sources surround the neutron star and contaminate mea- 
surements of the thermal emission from the point source. Also, sources with high 
E produce copious non-thermal X-rays from magnetospheric processes (see ii7.2|l . 
This emission, for example, completely overwhelms any thermal emission from the 
very young Crab pulsar, hence the availability of only an upper limit on its thermal 
emission, obtained off-pulse (Table [^J. Mitigation of both the above problems is 
best done by avoiding high E sources, even though they are often the youngest and 
hottest (see fT^ . 

An additional important difficulty is the degeneracy between the equivalent neutral 
hydrogen absorption toward the source, Nh, and the measured temperature. For 
an absorbed blackbody-like spectrum, the low-energy exponential cutoff due to Nh 
results in higher inferred T^ for smaller Nh, and vice versa. Generally it is necessar\ 



results m mgner interred 1^ tor smaller and vice versa, generally it is necessary 
to constrain T„°° vi a a contour plot in the T^/Nh plane (e.g. Halpern fc Wanel 1997t 
Zavlin et al.L 1X9991) . although frequently authors quote T^ for only a single assumed 
Nh (see Table El. 



7.3.3.2 Hot Spots and Cyclotron Resonance Scattering 
In many instances (see Table FT^ and ii7.5|l . the effective radius of the star as 
determined by either a blackbody or an atmosphere spectral fit is much smaller than 
the expected radius of a neutron star, regardless of assumed EOS. In at least one 
case, RX J1856.5— 3754, this has led to claims that the thermal spectrum th erefore 
implies the target is not a neutron star at all, but a quark star ( Drake et al.L l2002'). 

However, localized backheating of a neutron star surface by relativistic 
electron/positron pairs formed from conversion of 7-rays produced in the magneto- 
sphere in the polar-cap or outer-gap accelerators (see M7.2() can produce "hot spots" 
which emit thermal radiation. Such hot spots have areas much smaller than the stel- 
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lar surface and could in principle outshine surrounding thermal emission from initial 
cooling, thereby dominating the X-ray spectrum. This is the case for the thermal 
emission detected fr om mill i secoii d pulsars (e.g. PSRs J0030-I-0451, J0437— 4715, 
J2124— 3358; see e.g. iBeckeilEoOll for a review), which have ages of ~ 10^ yr hence 
should have cooled long ago. When hot spots dominate, studying thermal emission 
is mainly constraining the physics of the magnetosphcre, although upper limits on 
initial cooling emission could still be interesting. For example, the Vela pulsar is 
a strong 7-ray source, yet a very faint thermal X-ray source (see Tables 17.11 and 
I7.2|l : magnetospheric- acceleration-driven hot spots may well dominate the thermal 
spectrum, but its low effective temperature is still interesting (see Fig. 17.11)) . Thus, 
measured temperatures that are associated with small emitting areas are most con- 
servatively taken to be upper limits on cooling temperatures. 

Although one might expect stars with hot spots to have higher pulsed fractions 
than those without, this is not necessarily the case due to strong gravitational light 
bending expected in neutron stars (see H7.3.3.3(l . This and a less-than- favourable 
viewing angle can significantly reduce the pulsed fracti on seen from a s tar w ith hot 
spots, as has been argued for RX J1856.5— 3754 (e.g. Ransom et al. . 2002)) . The 
detection of a bow-shock Ha nebula around this source support s it b eing an off- 
beam rotation-powered pulsar fsee 17.51 van Kcrkwiik & Kulkarni| 2001). The same 
is true of other "isolated neutron stars" (see i)7.4() as well as of the "central compact 
o bjects" (see U7.5II 

R.udermai] l|2003l) has argued that cyclotron-resonance scattering of thermal pho- 
tons by electrons and positrons within a few stellar radii of an energetic rotation- 
powered pulsar results in a Planck-like X-ray spectrum that is quite changed from the 
seed thermal spectrum. A pair plasma is thought to surround the star, maintained by 
conversion of 7-rays from the star's polar-cap and/or outer-gap accelerators. This 
is because a thermal X-ray photon's energy, ^ 1 keV, is less than the cyclotron 
resonance energy at the stellar surface, eBh/27rmc = 11.6 keV for B = 10^^ G. 
Thus, as the photon passes through the magnetosphere, it goes through a resonance, 
and th e optical depth to cyclotron-resonance scattering becomes large. iR.udermanI 



l|2003l) describes the situation as if the X-rays are reflected back from the surface of 
a "lightly leaky Hohlraum container" and only escape after multiple reabsorptions 
and re-emissions from the stellar surface. Thus, he argues, measurements of ther- 
mal emission from neutron stars teach us only about the stellar magnetosphere, not 
about initial cooling. 



7.3.3.3 Pulsations and Polarization from Thermally Cooling Neutron Stars 
The pulsations from the surface of a cooling neutron star have low pulsed 
fraction and low harmonic content. This is in strong contrast to the high harmonic 
content and pulsed fractions seen from the magnetospheric component of the emis- 
sion, which is presumed to be produced well above the stellar surface (see t l7.2(l . 
Often young neutron stars exhibit emission from both mechanisms, as is clear from 
the energy dependence of their hght curves (e.g. Fig. l7.5|) . 

The low harmonic content results primarily from the presumably low variation of 
temperature over the neutron-star surface. However, pulsations are further reduced 
by the strong gravitational bending of light near the surface of the neutron star (e.g. 
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Pagelll995t IPsaltis et alil2000|) . This becomes important as the general relativistic 



compactness parameter p = Rc^ /2GM — > 1 (a Newtonian star is described by 
p oo). Due to bending, although the observer may be facing one side of the 
neutron star, she/he may still detect significant emission from the opposite side. 
Simulations of well defined infinite contrast spots on neutron stars show that the 
effect is so strong that eve n for jo = 4, the max imum pulsed fraction is less than 0.4 
for any size emitting spot ( Psaltis et all 2000() . 

Measured temperature is expected to be dependent on pulse phase, as the ther- 
mal conductivity of the neutro n-star surface depends on the m agnetic field strength. 



which varies over the surface ijPagei Il995t IPsaltis et all. l2000(l . However, generally, 



only phase-averaged temperatures are measured, both because of faint signals, and 
also because gravitational light bending makes it difficult to do otherwise. Clearly 
this complicates the interpretation since the averaged value depends on the distri- 
bution of temperature over the surface. However, this effect is likely to be small 
compared to other systematic effects. 

X-ray emission from the surface of a cooling neutron star endowed with a magnetic 
field at the surface of magnitude ^ 10^^ G is expected to be significantly polarized. 
This is because atmospheric opacities depend strongly on polarization when the pho- 
ton energies are much smaller than the electron cyclotron energy. The opacity to light 
with its electric field vector oriented perpendicular to the magnetic field is smaller 
by a factor of approximately the squared ratio of the photon energy to the cyclotron 
energy, if the rati o is sm all, relative to that for light with electric-field parallel to B 
I Lodenaual et al. Thus, radiation with electric-field vector perpendicular to 

B escapes from greater atmospheric depths where the temperature is higher, hence 
results in a higher flux. Po l arizations from localized regions on the surface should be 
very high, > 50% llKannoi Il97il : IPavlov fc Shibaiiovl , 1X97^ . However, an observer 
sees radiation averaged over a hemisphere, and even beyond due to gravitational light 
bending, whi ch is strongest for large M/ R. This will reduce the observed polarization 
fraction (e.g. Pavlov fc ZavlinL 20o3). Similarly, the different magnitudes and direc- 
tions of the lines of a dipolar B field at the n eutron-star surface might be thought 
to reduce polarization 1 Pavlov fc ZavlinL 2000() . However, the QED effect of vacuum 
polarization, in which the region surrounding the star is effectively birefringent, has 
recently been shown to counteract the GR and varying B-field effects, because the 
observed polarization direction of rays from the surface is correlated with the _B-field 
direction far from the surface; a ray bundle that we ob serve passes through only a 
small i?-field solid angle as it leaves the magnetosphere ( Hevl fc ShavivL 2002() . This 
results in very large net polarizations as a funct ion of pulse phase, and even phase- 
averaged polarization fractions of 10 — 20% ( Hevl et al. . 2003t Lai fc HoL BoOSi) . 
Measurement of polarization light curves should thus allow sensitive constraints on 
the B field, as well as on M/R and the magnetic geometry of the neutron star. 
Hence, thermally cooling neutron stars and magnetars (Chapter 14) are expected to 
be very interesting targets for future X-ray polarimetry experiments. 



7.3.3.4 Features in Neutron Star Thermal Spectra 

Features in the X-ray spectrum of neutron stars have been an astrophysical 
'Holy Grail" since their X-ray emission was discovered. This is because the presence 
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of lines of identifiable origin, with known inertial reference frame energies, would 
allow the determination of the stellar redshift, z, and hence, at least in principle, 
the mass-to-radius relationship via M/R — (c^/2G)[l — (1 + z)^^], which would 
constrain EOS. (The above is true only under the assumptio n that General Re l ativit y 
correctly describes the spacetime around the neutron star; Dedeo fc Psaltist 200.ll) . 
Atmosphere models with even as little as 0.05% metal abundance predict prolific line 
produc tion due to the many t ransitions available for the various atomic ionization 



stages i Raiaeopal et al. . 1997(l . Hydrogen atmosphere models, by contrast, predict 



few lines. The relevance of spectral features to constraining the neutron-star age is 
thus that (i) they can help determine the atmosphere composition which can help 
establish what continuum model to fit; and (ii) they could, in principle, distort the 
c ontinuum g iven poor spectral coverage or resolution. 

Sanwal et al. ( 2002_3), using Chandra's ACIS CCD detectors, observed two absorp- 
tion features, at 0.7 and 1.4 keV, in the X-ray spectrum of IE 1207.4-5209, in the 
supernova remnant PKS 1209-51/52 (G296 .5-H0.0). The nature of this object is 



not understood; see ^7.5\ Mereehetti et al.l |200 2al confirmed the existence of the 
lines using XMM-Newton. A subsequent long XMM-Newton exposure revealed an 
additional absorption feature at 2. 1 keV, as we ll as a possible feature at 2.8 keV (see 
Fig. mni Bignami ct al., 2003; De Luca et ail [20041. These suggest harmonics of a 
cyclotron line at 0.7 keV. A cyclotron fundamental at this energy implies a magnetic 
field of 8 X 1 0^° G for an electron line, o r 2 x 10^"^ G for a proton line. However, 
Sanwal et al. lj20 02b'l and Hailev & Moril l|2002[) argue against the cyclotron inter- 
pretation, and suggest possible spectral line identification for the most secure 0.7 
and 1.4 keV features. In any case, why such lines are not seen, in spite of many 
comparably sensitive observations of other neutron stars, remains a mystery. This 
source is discussed further in ii7.5l 

Very recently, features in spectra of three "isolated neutron stars" have been re- 
ported. These are discussed in ^7.4\ 



7.3.3.5 Neutron Star Ages 

The standard age estimator for rotation-powered pulsars assumes frequency 
evolution of the form (l — kft", where n is the "braking index," and fc is a constant 
that depends on the magnetic moment of the neutron star. Assuming k and n to be 
constant, the braking index can be determined from a measurement of the second 
time derivative of the frequency, Cl, via n = vCl/Cl^. The age derived by integrating 
the above differential equation is: 



{n - 1)P 



(7.4) 



where Pq is the spin period of the pulsar at the time it became a dipole rotator, 
generally presumed to coincide with the supernova event. For a simple vacuum 
dipole spin-down model, n — 3. For Pq <^ P and n = 3, Equation 17.41 reduces 
to Tc = P/2P, the often-used pulsar charact eristic age. For the fi ve pulsars for 
which a constant value of n has been measured llLvne et al lli988; Kas pj et al.Lll994l: 
Lvne et al l ll99fiHDeeter et all Il99i ICamilo et al.Ll200nl) . the observed values are in 
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Fig. 7.10. XMM-Newton spectrum of IE 1207.4-5209. The best-fit continuum 
curve is represented by tlie sum of two blackbody functions. Four absorption 
features appear at the harmonically spac ed energies of ~0.7 keV, ~1.4 keV, ~2.1 
keV and ~2.8 keV pe Luca et alll2004f l. 



the range 1.4-2.91. Though pulsars clearly do not rotate like perfect vacuum dipoles, 
the range of braking indexes is limited and observationally well constrained. 

The situation is less clear for the initial spin period, Pq. This can be determined 
from Eauation l7.4l if the age is known and n measured. This is only the case for the 
Crab pulsar, whose estimated Pq ~ 19 ms has led to the generally made assumption 
that Pq <^ P for all but the very fastest pulsars. However, the initial spin pe- 
riod distribution of neutron stars is not well predicted by theory, si nce the rotation 
rates of the cores of the massive progenitors are largely unknown lEndal fc Soial 
ll97S^ . Also, circumstances at core collapse could significantly affect the neutron- 
star spin independent o f the angular momentum properties of the progenitor (e.g. 
Spruit fc Phinnevl 1998() . As an exa mple, consider PSR J1811— 192 5 in the super- 



nova remnant Gil. 2-0. 3 (Fig. 171^1 iTorii et all Il999l: iKaspi et all boOlhl) . The 

supernova remnant age is well determined to be ~2 kyr, while Tc = 24 kyr, implying 
Pq — P 65 ms for reasonable n. In this case, Tc is certainly over an order of 
magnitude different from the true age. Si milar argumen t s hold for PSR J0205-)-6449 



in the historic supernova remnant 3C 58 I Murrav et al. . 2002() . 



7.3.4 Current Status of Observations Vs Neutron Star Cooling Theory 

Chandra and XMM-Newton have made by far the best measurements of 
neutron-star spectra thus far, and so offer the best opportunity for constraining 
models of neutron-star cooling. The current status of o bservations of thermall y 
cooling young neutron sta rs has recent ly been reviewed by Kaminker et al. ([2002^, 



Yakovlev et alJ l|2002l) and iTsuruta et al. (2002), and is summarized in Figure [T!H1 



There, the data points are actual temperature measurements or upper limits; those 
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Table 7.2. Observations of Thermally Cooling Young Rotation-Powered Pulsars 



7. 3 Thermal Emission from Cooling Neutron Stars 



31 




2 4 6 

Ig (t) (yrs) 



Fig. 7.11. Current best measured rotation-powered pulsar temperature as a func- 
tion of best age estimate (see Table 17.21 . "H" symbols denote that a B = 
hydrogen atmosphere model was used; otherwise, a _B = blackbody model was 
assumed. The cooling curves shown are for 1.35 Mq stars, as in Figure FTHl (solid) 
non-superfluid star with B — and no accreted mass (Fig. 17.81 left); (dot-dash) 
a combination proton and neutron superfluidity model (Fig. 17.81 middle); (dash) 
non-superfluid star with 5 = and accreted envelope mass 10"^ Mq (Fig. EH 
right). The curves are illustrative only and do not represent fits to the data. 



marked with an "H" are for hydrogen atmosphere models, while the rest are for 
simple blackbody models. The uncertainty from possible nebular contamination 
is not indicated; for sources with nebulae, the observations made with Chandra, 
given its superb spatial resolution, are superior to those with other telescopes. Age 
uncertainties are not indicated, but as discussed in t |7.3.3.5l can be an order of 
magnitude. For PSRs J0538-I-2817, B0656-I-14, and B0833-45 (Vela), we show both 
the blackbody and hydrogen atmosphere temperatures, to illustrate the substantial 
difference. The plotted data are summarized in Table [77^ 

A striking property of the data is its scatter. Although the above caveats must 
be applied, the plot remains log-log, and it is hard to explain all the data with iden- 
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tical neutron stars sha ring a common mass and hence EOS and interior structure. 
Kaminker et al. I (l2002l) argue that on the basis of the data together with models 
that account for most of the phenomena affecting cooling (they do not consider "ex- 
otic" particles, but do consider neutron and proton superfl uidity in detail), ther e 
are effectively three types of cooling neutron stars (see also iTsuruta et al . 20021) : 

(i) low-mass, slowly cooling stars where there is no direct Urea process occurring; 

(ii) medium-mass, moderately cooling stars in which some nucleon direct Urea pro- 
cesses occur, but for which proton superfluidity suppresses the bulk of the possible 
Urea processes; and (iii) high-mass, rapidly cooling stars for which suppression of 
direct nucleon Urea processes is weak. Of course the mass cuts corresponding to 
these different categories depend strongly on the EOS and on the poss ible types of 
superfluidity. As an example, for the EOS "A" of lVakovlev et alJ l|200llh . for neutron 
stars having age 4 kyr, for a 2p proton superfluid, the low- mass/medium- mass cut 
is 1.36 Mq, while the medium/high mass cut is 1.64 Mq. 

What is clear is that we have not yet seen any very rapidly cooling objects. How- 
ever this may be less of a statement about the interior physics as it is about the 
mass function of neutron stars, i.e. there may be very few high-mass neutr o n star s 



produced for purely astrophysical reasons. According to iKaminker et alJ l|200: 



the Vela pulsar (PSR B0833— 45, fit with a hydrogen atmosphere) and Geminga 
(PSR J0633+1746, fit with a blackbody), as weU as presumably PSR JG538+2817 
(with a hydrogen atmosphere) the coolest objects, are prime examples of moderately 
cooling stars for which some direct Urea reactions are occurring, and their masses 
could be determined assuming some EOS, a model of proton superfluidity (but see 
iTsuruta, et aTl . l2002l) and the dependence of the critical temperature for proton super- 
fluidity on density. Note that the hydrogen-atmosphere model fit to PSR B0656-I-14 
is generally ignored, as it yields a radius much greater than is possible for a neu- 
tron star. The apparently higher temperature objects, like PSR B1055— 52, must, 
in this case, be low-mass neutron stars, possibly affected by some internal heating 
mechanism or by polar-cap heating. The pulsar in 3C 58 recently received fanfare 
for its relatively low ter nperature in spite of having likely been born in 1181 AD, 
hence being very young ( Slane et al. . 2002|) . Though the Vela and Geminga pulsars 
nominally have lower repo rted temperatures hen ce might be argued as being closest 
to rapidly cooling objects 1 Yakovlev et al. . 20021) . in fact the 3C 58 upper limit was 
for a blackbody spectrum. A hydrogen-atmosphere fit to the data (as was done for 
Vela) could possibly bring the upper limit down to a very constraining level. Indeed, 
the Vela temperature when fit with a blackbody is well above that of the 3C 58 
pulsar. 



7.4 "Isolated Neutron Stars" 

A relatively recently discovered class of X-ray sources are believed to be 
thermally cooling neutron stars. Their nature is still not well understood. These are 
currently being referred to in the literature as "isolated neutron stars" (INSs). This at 
first glance may seem like a misnomer since this title does not distinguish them from 
isolated rotation-powered pulsars. Actually, as we describe below, it is possible that 
at least a subset of the INSs are conventional rotation-powered pulsars. However, the 
sources known as INSs are currently distinguished as showing no evidence for radio 
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pulsations, nebulosity, or accretion, and X-ray emission characterized by very soft 
spectra that are well described by a blackbody, with no apparent magnetospheric 
contributions. These very properties, as well as their proximity, make them enticing 
targets, since many of the observational difBculties inherent in using thermal emission 
to constrain the nuclear EOS that are described above in tl7.3l are not present. The 
name INS therefore does apply to rotation-powered pulsars as well, however the 
latter are a very special class of the INS population: long after rotation-powered 
pulsars have spun down past the "death line," beyond which they no longer produce 
observable magnetospheric emission (this typically takes ^--^30 Myr), "dead" neutron 
stars will still populate the Galaxy. INSs, if identified as this "dead" component, 
may well nu mber as high as 10 ^ in the Galaxy, much more than the ~ 10^ active 



radio pulsars l|Lvne et al 



Prior to the launch of the ROSAT satelli te, theoretical predictions following an 



idea first put forth by Ostriker et al. 1 197fl() were that several thousand otherwise 



dead neutron stars would be detectable in sof t X-rays by ROSAT. due to their accre- 
tion o f material from the interstellar medium l|Treves &: ColDiLllQQlllBlaes fc MadauL 
19931). In fact, wel l after the rise and fall of ROSAT, only 7 objects (Table O see 



Treves et al. . I2n0(l[ for a review) having the predicted properties (such as high X- 



ray to o ptical flux ratio) have been fo und in spite of very careful searches (e.g. 
Danneil Il998allj : Rutledge et al.L l2003|) . This is likely due mai nly to the much 



higher than expected average ne utron star space velocity (e.g. Lvne fc Lorimeil 



il9M: Hansen fc Ph inncv. 1997; Arzoumanian et al. . I2flfl2ll . which hampers stan 



dard Bondi-Hoyle accretion, although other explanations have been suggested (e.g. 



Pernaet alll2003l) . 



In fact, of the objects in Table [731 it now appears likely that none is accreting from 
the ISM; rather, they are visible because of thermal emission from the neutron-star 
surface. If so, these are unlikely to be extremely old neutron stars as they are too 
bright; rather they would have to be either off-beam conventional rotation-powered 
pulsars having ages ~ 10^ yr and magnetic fields 10^^ — 10^^ G, or even younger 
quiescent magnetars having magnetic fields 2-3 orders of magnitude higher. 

Although 4/7 of the INSs have now been shown to exhibit pulsations, none has had 
its spin-down rate measured. Their ages are thus not known, which hinders efforts 
to use them in cooling studies. However, in these sources, one hope is that spectral 
modeling, coupled with sensitive, high-resolution X-ray observations, could deter- 
mine their radii. If spectral lines were detectable as well, the stellar masses could 
in principle be known independently from the lines' gravitational redshift f ^7. 3. 3. 41) . 
thus providing a direct constraint on the neutron-star mass/radius relation and test- 
ing EOSs. Spectral modeling in these sources is greatly aided by their low equivalent 
neutral hydrogen column densities, typically 1 — 2 x lO^*' cm^^. Although these 
ambitious goals have not yet been realized, the recent dis coveries of absorption fea- 
tures in the spectra of three INSs is quite encouraging l|Haberl et al.L l2003t 
van Kerkwiik et all |^D4.) . 

The current situation is nicely illustrated by the example of RX J1856.5— 3754, 
which is the b est studied of the INSs. This object was discovered in ROSAT data 
I Walter et all Il996') and optic al observations with the H ubble Space Telescope re- 
vealed a faint blue counterpart ( Walter fc MatthewsL 1997(l . This suggested emission 
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Table 7.3. Reported Isolated Neutron Stars"" 



Name'' 


P 


ROSAT PSPC rate 


mv 


log fx/fopt 




(s) 


(cps) 






RX J0420.0-5022 


3.45 


0.11 




> 3.3 


RX J0720.4-3125 


8.39 


1.69 


26.8 


4.6 


RX J0806.4-4132 


11.4 


0.38 




> 3.4 


RX J1308.6+2127 


10.31 


0.29 


28.7 


5.0 


RX J1605.3+3249 




0.88 


27.1 


4.4 


RX J1856.5-3754 




3.64 


25.8 


4.4 


IRXS J214303.7+065419 




0.18 




> 3 



"Data from lTreves et al.1 (|200(t . iKaplan et alJ and lHaberJ (|2Qq1|2Qq3); 

see text for more information. 

''Source names ending with are uncertain members of the class. 



that is the Raleigh- Jeans tail of a kT ~ 50 — 60 eV blackbody. The HST detection 
enabled an astrometri c study of the source that detected both parallax and proper 
motion ijWalteil l200l|) . The source is evidently extremely close t o the Sun, having 
parallax 8.5 ± 0.9 mas I Kaplan et al. . 2002t Walter fc Lattimerl 2002|) or distance 
117 ± 12 pc. The proper motion of the source, if traced backward, appears to inter- 
sect the Upper Sco OB association, if it has a radial velocity of ^30 km s^^ I Walteri 
boOlilKaplan et aLlEooi I Walter fc Lattimerl Eool . Assuming the neutron star was 
born in the association implies an age of 0.4 Myr. 

Because of the proximity of the source and the chance for the discovery of spec- 
tral features, a 500 ks Director's Discretionary Time observation with the Chandra 
X-ray Observatory was performed in 2002. The X-ray emission was found to be un- 
pulsed, with stri ngent upper lirn i ts set on the pulsed fraction (<1.3% for frequencies 
less than 50 Hz; iBurwitz et al. The spectral results of the DDT observa- 

tion were a surprise. Not only were no features detected, but the resulting X-ray 
spectrum is extremely well modeled by a simple blackbody, with no evid ence for 
the atmospheric distortion so clearly predicted by theory (see H7.3.2.5I Drake et all 
2002t iBraie fc RomaniL l2002t iBurwitz et all l2003|) . In particular, the emitting ra- 
dius for a single blackbody fit is 4.3((i/117 pc) km, smaller than expected for any 
plausible neutron star. This led to speculation that the object is n ot a neutron star 
but rather a self-bound quark or strange star ijDrake et all 2002jl . A very recently 
revised parallax distance of 175 pc would imply a radius of 7.2 km, more in line 
with expectations for a neutron star (D. Kaplan, personal comm.). In any case, it 
is possible that the emitting region is simply a hot spot on a larger stellar surface; 
this is not completely precluded by the low pul sed fraction upper limits because of 



the strong effects of gravitational light bending l|Burwitz et all |2Q03) . Interestingly. 



the optical emission is too bright to be the low-energy extrapolation of the X-ray 
spectrum. A two-component blackbody model is necessary: the lower temperature, 
perhaps from the bulk of the surface and accounting for the optical emission, has 
kT°° < 33 eV, and the higher temperature, perhaps from a hot spot and accounting 



7-4 "Isolated Neutron Stars" 



35 



for the X-ray emission, has kT°° 63.5 eV . In this case, a strange st ar is not re- 
quired jPrake et all l2002t iBraie fc RomaniL l2002t iBurwitz et all l2003|) . Of course, 
the standard caveats regarding the use of blackbody models apply (see ii7.3.2.5|) . 
van Kerkwiik fc Kulkarnil l)200]|) discovered, using Ha imaging, a cometary shaped 



nebula having RX J1856.5— 3754 at its apex, similar to those detected around some 
conventional radio pulsars (see ii7.6|l . This suggests that RX J1856.5— 3754 is a 
conventional neutron star whose emission beam does not intersect our line of sight. 

Other sources in Table [731 are not yet as well characterized as RX J1856— 3754. 
We discuss each briefly, in order of the degree to which it has been s tudied. 

RX J0720.4-3125 is a 8.39-s soft- X-ray pulsar l|Haberl et al.Lll997(l plausibly iden- 
tified as an of f-beam conventional rotation-powered pulsar and not an accretor (see 



Kaplan et al.L 12003b . and references therein) . No spin-down rate has yet been de- 



termined, although int eresting upper l imits on P have been established, and imply 
that B < 3.5 X lO^^ G ijCronner et a,l .L l2004i^ . Its optical through X-ra.v spec trum is 
well modelled by two blackbodies (as in RX ,11856— 3754: lMotch et al. [ l2003l) plus a 
possible power-law component, similar to what is observed in young and middle-aged 
radio pulsars (e.g. PSR 60656+14). A large inferred velocity seeins to p reclude it 



being an accretor 1 Motch et al. . 2003|) . Intriguingly, Haberl et al. l|2004|) recently 



report the discovery of a phase-dependent absorption line it this pulsar's spectrum. 
The line has energy ~271 eV and an equivalent width that varies from —31 cV to 
—58 eV depending on pulse phase. The line appears best interpreted as cyclotron 
resonance of charged particles, presumably either electrons or protons. For the lat- 
ter, a 5 X 10^^ G magnetic field is implied. For electrons, the field must be 2000 
times smaller. 

RX J1308-I-2127 (RBS 1223) was identified bv lSchwope et all lll999l) in ROSAT 



data as a relatively bright INS candidate, based on it soft spectrum and the high X- 
ray to optical flux ratio made possible by the good posi tional localizat i on, lilaberl et all 



2003() showed the source is pulsed with P = 10.3 s. iKaolan et al.l l|2002fl reported 



a possible optical counterpart having flux in excess of the extrapolation of the X- 
ray blackbody spectrum, as in RX J1856-3754 and RX J0720.4-3125. A broad 
absorption line feature has been reported in this source JHabe rl et al.[ 20031). and 
interpreted as a proton cyclotron line in a > 10^^ G field. Hab erl et al.l l|2003l) re- 
ported the discovery of a broad absorption feature in this object's X-ray spectrum. 
The line energy is at ~300 eV and has equivalent width ~ —150 eV. This line is 
interpreted by the authors as being due to cyclotron resonance absorption, in the 
lO^'^ — 10^** G range if protons, and in the lO^'' — 10 ^^ G range if electro ns. 

RX J1605.3+3249 (RBS 1556) was identified by iMotch etaD l|l999|) as being a 
possible INS, because of the soft source spectrum (kT ~ 9 2 eV), the lack of X- 



Kaolan et al] l|2003ah 



ray variability, and the absence of any optical counterpart, 
reported the detection of a faint, blue optical counterpart, which they argue confirms 
the INS interpretation. The optical counterpart has flux that is over an order of 
magnitude above the Rayleigh- Jeans extrapolation of the X-ray blackbody. This is 
like, though somewhat larger than, the excesses seen in the other sources, and argues 
either for the X-rays being from a hot polar cap, or, for the optical emission having 
a non-thermal component. Either is consistent with an off-beam rotation-powered 
pulsar. Very recently, Ivan Kerkwiik et alJ l|2004j) report the detection of a broad 
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absorption feature at 450 eV. As for RX J0 720.4-3125 and RX J1308+2127, the 
origin of the hne is as yet unclear, although Ho fc Lai lj2004|) show on the basis of 
approximate calculations that all can be plausibly explained by proton cyclotron 
absorption. 

as an interesting INS can- 



RX J0806.4-4132 was reported bv lHaberl et al 



didate. It was seen in two RO SAT AW Sky Survey observations at the same flux, and 
having a soft spectrum, with kT ~ 78 eV. Haberl et al. 1 19981) detected no optica l 
counterpart, indicating a high X-ray to optical flux ratio. ^Haljerl fc ZavlinI l)2002l) 
using XMM-Newton observations of this source, report a possible 11.4 s periodicity 
at the 3.5cr level having a small ^^6% pulsed fraction, and a blackbody temperature 
of kT = 94 eV. The improved XMM-Newton position should allow deeper searches 
for a possible optical counterpart. 

RX J0420. 0—5022 is a iZ05'j4T' soft-spectrum source showing evidence for a 22.7 s 
periodicity, a high X-ray to optical flux ratio, no apparent variabili tv, and a soft 
spect rum that is well described by a blackbody of temperature ^^57 eV (jHaberl et al. 



19991) All argue for it being an INS, although the period is much longer than has 
been seen for any rotation-powered pulsar. 



7.5 "Central Compact Objects" 

There are several putative compact objects that have been studied exten- 
sively and which have potential to constrain neutron-star cooling theory fi l7.3|l . yet 
which are poorly understood. By definition, these sources, known as "Central Com- 
pact Objects" (CCOs), are found near the centers of supernova remnants (SNRs), 
although this could well be a selection effect resulting from SNR X-ray studies. The 
CCOs share the following distinguishing properties: unusual X-ray spectra that, flt 
with blackbodies, and, in some cases, hydrogen atmospheres, imply very small radii, 
and very high effective temperatures (see Table 17^ . The small radii are inconsis- 
tent with these objects being neutron stars for any EOS, while the temperatures 
are inconsistent with any neutron-star cooling models. Other distinguishing char- 
acteristics are high X-ray to optical luminosity ratios, no evidence for pulsations 
(except in one case ~ see below), no evidence for a wind as seen in conventional 
young rotation-powered pulsars fi i7.6(l . and no evidence for any companion star that 
could be powering the X-ray emission via accretion. We discuss th e five be st studied 
CCOs briefly below. For a more detailed review, see Pavlov et alJ 1 2002a ). 

The "first light" of Chandra targetted the 300-yr old supernova remnant Cas- 
saeopia A, and revealed, for the first time, the presence of an X-ray point source 
near the center fFig. 17111 iTananbaumL [1999 *). This object, CXO J232327.9-K584843, 
seems likely to be the long-sought-after compact object formed in this oxygen-rich 
remnant, clearly a result of the supernova explosion of a massive star. The spectral 
properties of this source, as well as the absence of any detectable pulsations, make 
it ver y different from canonica l young neutron stars like the Crab pulsar. Specif- 
ically, Pavlov et al. (200(1) and Chakrabartv et al. 1 2001 ) showed that for either a 
simple blackbody or hydrogen atmosphere model, the temperature of the emission 
{kT = 0.25 — 0.35 keV) is very high, and the emitting area {R = 0.2 — 0.5 km) too 
small to be consistent with surface thermal emission from initial cooling. A power- 
law fit yields a photon index (F 2.6 — 4.1), significantly higher than those seen in 
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Table 7.4. Central Compact Objects''^'' 



Name 


SNR 


P 




poo 
^bb 












(keV) 


(km) 


(keV) 


(km) 


CXO J082157.5-430017 


Pup A 




0.4 


1.4 


0.2 


10 


CXO J085201.4-461753 


G266.1-1.2 




0.4 


0.3 


0.3 


1.5 


RX J121000.8-522625 


G296.5+10.0 


0.424 s 


0.25 


1.6 


0.1 


1 


CXO J161736.3-510225 


RCW 103 


6 hr ? 


0.4-0.6 


0.2-1.6 


0.3 


1-8 


CXO J232327.9+584843 


Cas A 




0.5 


0.5 


0.3 


1 



"See lPavlov et alJ ^2002a^ ■ IPavlov et alJ ll2002lj) and references therein. 
''Note: kT^ and R'^ refer to blackbody models, while kT°° and 7?°° refer to 
magnetic hydrogen atmosphere models. 



pulsars. The difference is further evidenced by the absence of any synchrotron nebula 
around the object (F ig. 17.121) . as well as the absence of any radio emission from the 
source I McLaughlin et al., 20^. The very high X-ray to optical flux ratio, however, 
does implicate a compact obiect and rende r an accreting binary scenario problem- 
atic IjKaplan et al. The source spectrum is consistent with 
it being an "anomalous X-ray pulsar" (see Chapter 14), that is, endowed with an 
ultrahigh magnetic field (10^"' — 10^^ G), or with it being a conventional B ^ 
neutron star having a coolin g iron surface an d hot hydrogen or helium polar caps 
I Pavlov et al . 2000: Ch akrabartv et al.U200lll. More recent XMM-New ton observa- 
tions of the source support these conclusions I Mereghetti et all 2002bj) . However, if 
it is an AXP, why it should be so faint and show no pulsations, and, if it has hot 
polar caps, why they should be so hot in the absence of any evidence for the pulsar 
mechanism, is a mystery. 

CXO J082157.5-43001 7 is located in the Puppis A SNR, thou gh signifjcantly off- 
center JPetre et allllQsl . From ROSATand ASCA observations lPetre et allll99(il: 



Zavlin et al.L Il99^ . the object's X-ray spectrum is well modeled by a blackbody, 
however the best-fit radius, 1.4 km, is much smaller than that of neutron star. A 
magnetized {B ^ 6 x 10^^ G) hydrogen atmosphere model yields a more neutron- 
star-like radius o f 10 km, and a tem perature more in line with that expected from 
cooling emission (|Zavlin et al.lll999^ . However, neither pulsations nor evidence fo r 
a surrounding nebula has been seen (|Gaensler et al. . 2000al: Pavlov et al. . 2flfl2al) . 
arguing that if this is a rotation-powered pulsar, it must have very low E. This 
argues against the heated-polar-cap interpretation, as in Cas A. 

SNR G266. 1 — 1.2, also known aff ectionately a s "Vela Jr," is in the direction of the 
southeast corner of the Vela SNR l|AschenbachI Il998l) . The possible detection of a 
radioactive ^"^Ti line suggested an age of ~700 yr, and a distance of ^200 pc. How- 
ever, subsequent ASCA observations yielded a Nh more consist e nt wi th a distance 



of ^ 1 kpc, and an age of several thousand years l|Slane et all l200l|) . The point 



source CXO J085201. 4— 461753 is near the center of the SNR, and h as an X-ray to 



optical flux ratio only consistent with a compact object identification (jPavlov et al 
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Fig. 7 .12. "First light" Chandra X-ray image of the supernova remnant Cas A 
(after l^nanbaum, 1999). The image is 6' on a side. The previously unknown 
point source at the center of the remnant is obvious. 



20011) . This source is no t seen to pulsate and sh ows no evidence of a pulsar wind 



nebula in Chandra data I Kargaltsev et al. . 2002|) . The source spectrum, fit with a 



blackbody model, yiel ds a radius of only ^ 0.3 km. Atmosphere models do not yield 



20021) . 



more reasonable radii ijKargaltsev et al. 

RX J121000. 8— 522625 is a compact source near the center of the large SNR 
G296. 5+10.0 (also known as PK S 1209-51/52). The sour ce was first noticed i n Em - 
stein IPC observations made bv lHelfand fc Beckeil ^98^. iMereehetti et alJ l|l996|) 
used ROSAT to show that a blackbody model fit the spectral data well. They fur- 
ther demonstrated the absence of any extended emission, and set an optical limit of 
y ^ 25. They did not detect any radio emission from the so urce at 4.8 GHz, and 
Kaspi et al.l 1*1996') saw no pulsed radio emission at 436 MHz. Vasisht et al. 



and lZavlin et al.. (^199811 presented ASCA observations of the source which again sup- 
ported the neutron-star interpret ation. Chandra obse rv ations suggested a p eriodicity 



in the source, with P — 424 ms jZavlin et al.U2000(l . IPavlov et alJ lj2002j) . using a 



second Chandra observation, reported a spin-down age of 200-900 kyr, much larger 
than the 3-20 kyr age inferred for the SNR. However, subsequent observations have 
not found simple spin down as in rotation-powered pulsars, unless the source is an 
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active glitcher l|Zavliii et al. , 2004|) . This is currently not understood. The source 



spectrum is thermal , and the best-fit co ntinuum blackbody model yields a small 1.6- 
km emitting radius 1 Pavlov et al-lbOO^ . This source also shows absorption features 



in its X-ray spectrum fFig. mUl ti7.3.3.4|l . 

IE 161348—5055, the central source i n the supernova remnant RCW 10 3 was 
the fir st proposed cooling neutron star ijTuohv & Garmirc, 1980). Gotthclf et al.l 
ljl999|) showed that the source is X-ray variable by as much as an order of magnitude 
on a time scale of years. A ~6 hr p eriodicity was t entatively reported in data 



from Chandra and ASCA observations l|Garmire et all 12000.) . and a 2002 Chandra 



observation sho wed a clear 6.4-hr per iodicity and apparently phase correlated X-ray 
flux variations feanwal et al. . 2002a|) . This strongly suggests that this object is, in 



fact, an accreting neutron star having a very low-mass companion. 

7.6 Pulsar Wind Nebulae 

7.6.1 Physical Overview 

Only a small fraction ( < 10%) of the spin-down energy of a young pulsar is 
converted into observable pulsed emission ( WI .2\ . It is generally accepted that most 
of the energy leaves the pulsar's magnetosphere in the form of a magnetized wind 
I Michel 19691) . In the ideal (E ■ B = 0) MHD approximation to the aligned rotating 



magnetic dipole jGoldreich fc ,Tulianlll965lContoDoulos et al.L ll999>). it appears that 
the poloidal field (Bp, the field perpendicular to the toroidal component B^) goes 
from a nearly dipolar structure within the light cylinder to a split monopole struc- 
ture outside the light cylinder (i.e. Bp = B^ (x. r~^, a purely radial field whose sign 
abruptly changes at the equator). The toroidal component i?^, which is small out 
to the light cylinder grows rapidly outside the light cylinder so that for radii greater 
than a few times Rlc, >> Bp. "Cold" (i.e. non-radiating) charged particles 
flow outward with this magnetic field forming a magnetized wind which is ultimately 
accelerated to very high energies (7 >> 1). This highly relativistic magnetized wind 
eventually interacts with the surrounding medium, and emits synchrotron radiation 
from radio to 7-ray wavelengths. This synchrotron emission caused by the pulsar 
wind is what is generally meant by the term pulsar wind nebula (PWN). However, 
the winds from rapidly moving pulsars may produce a bow shock in the ISM caus- 
ing Ha emission in what is essentially a thermal pr ocess. Optical fil aments in the 



Crab nebula also appear to be thermal emitters fee. iFesen fc Kirshner . 1982). and 
in principal there could also be thermal X-ray emission at the outer edges of young 
PWNe, although this is difficult to distinguish from the surrounding supernova rem- 
nant (SNR) and has yet to be observed unambiguously . 

The details of the structure and luminosity of the PWN should depend on the 
pulsar's spin-down energy history and space velocity as well as the density profile of 
the surrounding medium. It may also depend upon the magnetic inclination angle 
of the pulsar, although how the wind properties depend on this angle is poorly 
understood and it is possible that the deviation from the aligned rotator model is 
negligible in terms of overall growth and energetics of the PWN. What we observe 
is dependent on the observer's viewing angle, with Doppler boosting being a non- 
negligible factor in the surface brightness distribution of the PWN. 
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Since a pulsar is born in a supernova, the initial environment that the pulsar wind 
encounters is the (nearly) freely expanding supernova ejecta. The pulsar may be 
born with a large space velocity and overtake the expanding ejecta shell within a 
few tens of thousands of years, after which the wind environment is probably the 
ISM There appears to be a minimum spin-down energy necessary to create bright 
PWNe with a s ignificant drop in PW N emission efBciency from pulsars wit h log 
i;<36 jPrail fc Scharringhausenlll997HGaensler et allbOOObl: ICotthelA l2003j) . and 



generally only the youngest isolated pulsars are observed to have PWNe. However , 
millisecond pulsars can also have a substantial E, and both Ha I Bell et al. . 1995(l 

and X-ray ( Stappers et al. , 2003(l nebulae have been observed aro und them. 

The standard theoretica l picture for PWNe from young pulsars llPacini &: Salvati 
197^ iRees fc CrninL IT974I Eevnolds Che^^ Il98i iKenne] fc Coroniti Il984ah 
is of a synchrotron bubble being blown at the center of an expanding supernova 
remnant. Initially, the outer edge expands supersonically into the supernova ejecta, 
so is presumably bounded by a forward shock {Rp of Fig. 17.131 . The corresponding 
reverse shock (Rt) is near the center of the nebula, where the cold relativistic wind 
from the pulsar is terminated. (Note that the wind is cold only in the sense that 
kT << (7 — l)mc^ and hence in the co-moving frame of the bulk flow the particles 
have very little time to radiate before encountering the shock.) The bulk flow en- 
ergy is then converted into random particle motion with a power-law distribution of 
particle energies. The wind continuously injects high-energy electrons and positrons 
as well as magnetic field into the bubble. The particles gyrate in the magnetic 
field, emitting synchrotron radiation from the radio through soft 7-ray regions of the 
spectrum. 

7.6.2 Observational Properties of PWNe 

Except for a few cases (notably the Crab and Vela nebulae), we have data 
on PWN emission only at centimeter wavelengths an d /or in X-rays be l ow ^ 10 keV. 
The radio properties of PWNe have been outlined by Weiler fc Shaved ( 1978|) among 



others. A source can be identified as a radio PWN if it has an amorphous (i.e. non- 
shell) morphology, flat, non-thermal spectrum (energy spectral index a ~ 0.0 — 0.3, 
where flux S — i/^"*), and high ( > 10%) fractional linear polarization. The X-ray 
emission also has a spectrum described by a power law, except the photon index 
r = 1 -I- a is generally between ~ 1.5 — 2.2. Hence there is a need for one or more 
spectral breaks such that the net change in spectral index is Aa = 0.5 — 0.9. The 
X-ray nebula tends to be smaller than the radio nebula, and generally has a more 
defined morphology, sometimes taking the form of a thick torus with perpendicular 
jets. 

The efficiency with which the spin-down energy is converted into synchrotron ra- 
diation can be quite high. However, it is often difficult to determine since the ob- 
servations tend to be restricted to two narrow bands. In general, one can find four 
observed quantities, with associated errors, quoted in the literature (or minor vari- 
ants): radio flux density at ~ 1 GHz, S'ighz, radio energy spectral index a^, 1-10 
keV X-ray flux (corrected for interstellar absorption) , and X-ray photon spectral 

* In radio work, the spectral index is often defined as 5" = i^" . However we use the minus sign here 
to be consistent with the X-ray definition. 
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Fig. 7.13. Density profile of a spherically symmetric PWN expanding within a 
pre-Sedov SNR showing the PWN termination shock Rt, the PWN forward shock 
Rp, the SNR shell reverse shock Rn, the contact discontinuity between the SNR 
eiecta and the ISM Rap, an d the SNR shell forward shock Rf- Adapted from 
ivan der Swaluw at all l2001f) 



index F = 1 + c^a;- From these, one can infer a break frequency where the spectral 
slope changes: 

"'-(IS-) ■ (") 

where Sx and Sr are the flux densities at Vx and Vr, respectively. Possible physical 
meanings of this break frequency will be discussed in i l7.6.4l At centimeter wave- 
lengths, the spectral index ctr < 1, hence the bulk of the energy emitted in this 
band is at higher frequencies. The radio luminosity Lj. therefore depends on what is 
chosen as the upper frequency limit v^i of the radio band: 

i, = i.ixl03O^i^52^^^^ergs-i, (7.6) 
1 — ar 

where d is in kpc, 5ighz is in Jansky, and Vui is in GHz. A natural upper limit would 
be the break frequency, but that is usually not well known due to uncertainties in the 
radio and X-ray spectral indices and the large extrapolation. Often, an upper limit of 



100 GH z is randomly chosen in order to compare different nebulae (jGaensler et al, 
2000bl) . 



The X-ray luminosity Lx can be defined as the total power emitted from the break 
energy up to a high-energy cutoff Ec ■ 



Lx = 1.1 X 10-:^^^g£;-^^ergs-, (7.7) 
where Fxi2 is the 1-10 keV flux in units of lO^^^ergcm^^ s~^ and Ec and Eb = hi^b 
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are in keV. If the photon index F < 2, the X-ray luminosity depends critically on 
the high-energy cutoff, which is only known in a very few cases. For the Crab and 
Vela PWNe, it is in the few to tens of MeV range. If F > 2, then the bulk of the 
energy is emitted below the X-ray band, peaking near the break frequency. Note 
that the X-ray luminosity is usually much higher than the radio luminosity. Since 
the synchrotron lifetime of the X-ray emitting particles is often much shorter than 
the age of the nebula, the properties of the X-ray emission can be highly dependent 
on the current spin-down energy of the pulsar. The radio emitting particles generally 
have synchrotron lifetimes longer than the age of the nebula, and are therefore an 
indication of the total number of lower-energy relativistic particles injected into the 
nebula throughout its lifetime. 

All known pulsars with detected X-ray or radio PWNe are listed in Table 17.51 
Table Ol lists other X -ray or radio PWNe from which pulsations have not been 
detected, while Table lTTzl lists pulsars with Ha PWNe. Note that the isolated neutron 
star RX J1856.5— 3754 also has a Ha nebula associated with it (see >I7.5|I . A source 
is considered a PWN if it has one or more of the following observational signatures: 
a known pulsar embedded in a radio and/or X-ray nebula which appears to be 
morphologically related to the pulsar; an isolated, fiUed-center radio nebula with 
the properties mentioned above; a composite supernova remnant defined as a SNR 
shell with cither a central radio component whose spectrum is flatter than the shell, a 
compact but extended central X-ray component with a power-law spectrum, or both; 
an X-ray point source with a very high ratio of X-ray to optical flux and associated 
extended X-ray emission with a power-law spectrum. The X-ray morphologies are 
quite varied, but can be loosely placed within the following categories: torus -I- weak 
jet (T), jet /trail dominated (J), and mixed or uncertain (M). In the mixed case, 
sometimes there appears to be a bow-shock or torus as well as a bright trail or jet- 
like component, but sometimes there is just an indistinct blob. Although a dominant 
toroid is generally assumed to be the typical structure, only a few of the brightest 
sources have been clearly determined to have toroidal morphologies. A reference for 
X-ray morphology is given where available, otherwise for the radio morphology. 

There seems to be a change in PWN morphology from toroidal to more of a jet 
or trail when the pulsar spin-down energy drops below logE ~ 36.5 (in erg s""'^). 
This probably simply reflects the on-average greater age of these pulsars so they 
have moved substantially from their birth sites. The radio and X-ray PWNe largely 
disappear once the spin-down energy drops below logi? ~ 35.5. Those few sources 
listed as being PWNe whose pulsars' spin-down energies are below this tend to have 
barely detectable X-ray PWN. The fact that they are observable at all, along with 
the existence of Ha nebula seen around pulsars whose \ogE < 34.0 attest to the 
persistance of the relativistic wind generation mechanism long after the effects on 
the surrounding medium are readily apparent. 

7.6.5 Evolution of a PWN in a SNR 

The standard evolutionary model of PWNe is a spherical bubble being blown 
at the center of a spherical SNR shell (Fig. I7.13|l . Observationally, this picture 
does not accurately represent the true state of affairs. Many PWNe, including the 
"canonical" example of the Crab, do not have observable SNR shells, and so there 
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"X-ray Morphology. T = torus -|- weak jet, J = jet or trail dominated, M = 

mixed or uncertain. A ? implies that no X-ray nebula has been clearly detected 
yet, but if there is a clear radio or Ha morphology, that is indicated in 
parentheses. 

'indicates if there is an associated SNR shell. O means it is outside the shell, ? 
means there is a nearby structure that could be interpreted as a related SNR 
shell. 

'^Indicates coincident 7-ray source. ? implies outside nominal 7-ray error box but 
possibly still related. 



is no observational evidence for the expanding blast wave caused by the supernova. 
However, presumably there is an external blast wave although the surface brightness 
of the shell might be very low. PWNc appear to be highly asphcrical, sometimes with 
thick torii, sometimes dominated by narrow, jet-like features. How this asphericity 
may affect the evolution of the PWN is not well understood, although it is plausible 
that the mere addition of a filling factor to the standard equations could adequately 
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Table 7.6. PWNe With No Detected Pulsar 
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See footnotes to Table 1731 



account for the morphological variations ( Roberts et all [2003). Pulsars are generally 
high- velocity objects, so they rapidly leave their birth site, eventually overtaking the 
supernova blast wave. The interaction of the pulsar wind with the swept-up shell 
material can reenergize the shell, which seems to be occuring in the case of CTB80 
I Shull et al. . 1989t van der Swaluw et all 2003aj) . To reach the shell usually takes 
a few tens of thousand years, and so for the first few thousand years assuming the 
pulsar is near the center of the SNR is reasonable. Well before the pulsar reaches the 
SNR shell, it will overtake the outer edge of its original bubble, severely distorting 
the shape, changing from a spherical/toroidal structure to more of a bow-shock. 
Since t he forward shock of the FWN is generally accelerating within the expanding 
ejecta I Reynolds fc Chevalieil 1984|) . this will usually not happen until after the 
PWN forward shock encounters the SNR reverse shock. Analytical and numerical 
models of t his tra nsition suggest it occurs at roughly half the pulsar crossing time of 
the nebula l|van der Swaluw et al., 2003b). 

Therefore, for the first few thousand years of a pulsar's life, the assumptions (other 
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Name 


Pulsar 


log(i5) 
(ergs 1) 


d 

(kpc) 


V 

(kms-i) 


Reference 
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1.9 


260 


Jones et al. 02 
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B1957+20 
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94 
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3.2 
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than sphericity) of the standard mod els are probably valid. The evolutionary picture 
has b een described analyti cally (e.g. IPacini fc Salvati'. 1973; Reynol ds fc Chevalie 
Il984|) and numerically (e.g. van der Swaluw et al.. . .2001: ,.Blondin et alTT20nih . Since 
the content and structure of the pulsar wind are poorly understood, it is generally 
assumed that a spherical outflow of energy equal to the spin-down energy of the pul- 
sar (the observed pulsed electromagnetic radiation being a negligibly small fraction 
of the energy) continuously flows outwards from the pulsar magnetosphere. This 
energy is split between the kinetic energy of the bulk particle flow and the magnetic- 
field energy. At some radius which is small compared to the total radius of the 
nebula, the electron/positron component of the bulk flow is randomized and takes 
on a power-law distribution of energies. If this transition is caused by a standing 
shock, then for efficient shock acceleration of the particles to occ ur, the flow energet- 
ics mus t be dominated by the par ticle componen t at th e shock. Kennel fc Coroni'til 
" l984a|^ . expanding on the work of IPees fc GunnI 1 19741) . constructed a steady state 



MHD model of the Crab nebula, assuming a positronic wind terminated by an MHD 
shock whose downstream properties are determined by boundary conditions imposed 
by the size. X-ray luminosity and expansion velocity of the nebula. The upstream 
flow properties can then be determined from the Rankine-Hugoniot relations. The 
ratio of magnetic to part icle energy flux a in the unshocked wind of the Crab is 
inferred to be cr ~ 0.003 (|Kennel fc CoronitiL ir984b). This leads to the well-known 
"sigma problem" of pulsar wind theory, since mechanisms for accelerating the parti- 
cles within the magnetosphere generally require a > 1 (see i l7.2|l . This suggests that 
the bulk of the acceleration occurs somewhere in the wind zone outsi de of the pulsar 
magnetosphere, for which there is no generally accepted model (e.g. IMelatosl 1998t 
Kirk fc SkiajraasenllioO^ . 



Although shock theory requires the ratio of magnetic to particle energy just down- 
stream of the wind termination shock to be small, it is expected that the magnetic- 
field will grow rapidly towards equipartition in the downstream flow. This allows 
the assumption that the magnetic and particle energy densities are roughly equal 
throughout the bulk of the PWN. For the flrst few days after the supernova explo- 
sion, the magnetic field density rapidly increases within the nebula, but then the 
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expansion of the nebula causes the magnetic-field strength to decrease with time 
despite continuous injection of magnetic energy. 

The subsequent evolution of the PWNe can then be separa ted into two eras of the 
pulsa r spin-down E and two eras of the SNR expansion (e.g. Reynolds fc Chevalieil 



11984 - For times t < t = Po/2Po, the initial spin-down time of the pulsar (note 
that this is probably the case for systems where the true age t << Tc = P/^P, the 
observed characteristic age of the pulsar), energy will be injected at a nearly constant 
rate {E ~ Eq), and the radius of the PWN expands as i? oc t^^^ while the average 
magnetic field strength will decrease as -B cx t^'^. For large initial spin periods 
(Pq > 50 ms) , the time of nearly constant energy input can last for several thousand 
years. During this time, the outer edge of the PWN bubble is moving supersonically 
compared to the surrounding medium, and so a forward shock should form, heating 
the swept-up ejecta which is in a thin shell around the PWN. The total mass of 
the swept-up material is initially fairly small, and so the thermal luminosity of this 
material may not be very large. At later times > r or t ~ Tc) when E << Eq, the 
expansion slows to i? (x i if the SNR is still in its free expansion phase. 

The blast wave of the surrounding SNR shell initially expands almost freely into 
the ISM (or into a surrounding bubble blown by the wind of the pre-supernova 
progenitor starjchcvalicr, 1982), shock heating the swept-up material. As the swept- 
up mass builds, a reverse shock forms and begins to propagate towards the center 
of the SNR, heating the ejecta. This reverse shock will reach the center of the SNR 
when the swept-up ISM mass is much greater than the e jecta mass, m arking the 
beginning of the so-called Sedov- Taylor phase of the SNR ijSedovl 19591) . Until this 
time, there are four shocks in the system: the forward shock of the SNR blast wave 
propagating into the ISM, the reverse shock of the SNR blast wave propagating back 
into the freely expanding ejecta interior to the shock, the forward shock of the PWN 
bubble, and the PWN reverse shock/wind termination shock fFigure lV.lSII . When 
the SNR reverse shock encounters the PWN forward shock, numerical simulations 
show the PWN is initially compressed, and then its radius oscillates a few times 
over the next several thous and years before settling dow n into a more compact and 
probably distorted shape ( van der Swahiw et al.L 12001 ) . During the compression. 



the magnetic field is enhanced causing an increase in the s ynchrotron luminosit y 
and filaments can develop due to Rayleigh- Taylor instabilities l|Blondin et al.ii200ll) . 
If the pulsar wind luminosity is still near its initial luminosity (E ~ Eq) after this 
compression, the PWN begins to grow again, but now the expansion is subsonic with 
R oc /:ii/i 5 jf^ other hand. E < < Eq then the expansion of the PWN goes as 

R oc t^/^ llRevnolds fc Chevalieil. Il98l . 



7.6.4 Emission Spectra of PWNe 

The relativistic particle gas that fills the PWN bubble loses energy through 
synchrotron radiation and adiabatic expansion. Since the observed radiation spec- 
trum has the form of a broken power law, a power-law electron energy spectrum is 
suggested, with the particle spectral index s (defined as N{E) cx i?"") related to the 
radiation spectral index a by s = 2a -I- 1. The simplest model is of a single electron 
population injected into a unifor mly magnetized n ebula with a single spectral index 
S and a high energy cut-off £;ec l|Chevalieil [2000(1 . Synchrotron radiation is usually 
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the most important cooling process for the X-ray emitting electrons, with the typical 
timescale tc for cooling dependent on the electron energy and the magnetic field. At 
high energies, « t, the age of the PWN, and an approximate steady state can 
be assumed between the injection and cooling of the electrons. In this regime, the 
average particle spectrum in the nebula is 5' + 1 which corresponds to a steepening 
of the observed spectrum by Aa — 0.5. A break in the particle spectrum occurs at 
an energy where tc ^ t, resulting in a break in the observed photon spectrum at: 

where B is in Gauss and t is in years. Note that since B decreases with time, this 
break frequency actually increases as the PWN ages. At low energies [tc >> t), cool- 
ing losses are negligible, and the particle spectral index is simply that of the injection 
spectrum. At late times when E << Eq, the initial cooling of the large number of 
particles injected at early times can leave an imprint on the particle spectrum in 
the form of a separate break below j/f, whose frequency decreases with time due to 
losses from the adiabatic expansion of the nebula. This res ults in an intermediate 
parti cle spectral index S* > s > 5 -I- 1 between the two breaks ( Revnolds fc Chevalieil 



19841) . The passage of the reverse shock causes a brightening of the nebula and an 



enhancement of the magnetic field. The high-energy electrons all rapidly lose their 
energy, resulting in a spectral cutoff of the electron population. After the reverse 
shock passes, the low-energy spectrum is dominated by the old electrons with en- 
hanced emission, while the high-energy spectrum is from newly injected particles, 
with a middle range below the post-reverse-shock cooling break which could have a 
fairly steep spectrum, i.e. Sm > S + 1. 

It is tempting to identify the radio emission with the uncooled s = S spectrum 
and the X-ray emission with the cooled s = S + 1 spect rum. However , this predicts 
Aa = 0.5 and in many cases it is significantly greater i Woltier et al. , 1997j) . If the 



X-ray or radio emitting particles were in the adiabatically cooled regime, Aa < 
0.5. The only time a single injection spectrum would fit this situation is after the 
reverse shock passage if the cooling break is above the soft X-ray band. In this 
case, there should be a spectral hardening somewhere in the hard X-ray or soft 
7-ray band. Another possibility is that the injection spectrum has changed over 
time. However, a correlation between PWNe X-ray spectral indices and E showing a 
softening with increasing E has been claimed ICotthcJi, 2003), implying the spectral 
index is decreasing with age. This would again lead to Aa < 0.5. A third possibility 
is that more than one electron population is being inje cted into the n ebula, one 



dominant at low energies and a second at high energies ijBandiera et al.l i2002,) . It 
is now clear that many, if not most, PWNe consist of both equatorial outflows and 
polar outflows, the latter in the form of coUimated jets. It might be expected that the 
two flows would inject particles into the PWN with different spectral slopes. Such a 
situation would often require the polar flows to be relatively underluminous in the X- 
ray region. A more funda mental problem is the flatness of the PWN spectra in radio. 
Both theoretical models I Kirk et al. , 2000(l and observations of other astrophysical 



shocks which produce radio emitting electrons (eg. SNR shells and AGN jets) tend 
to produce emission spectral indices of a ~ 0.6 and then cool to a ^ 1.1. The 
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only workable mechanism for producing the radio spectra of PWN currently in the 
literature is one developed for the Cra b where a low- energy cut-off which is rapidly 
cooled during the early spin-down era I AtovanL 1999^ . However, this requires some 
fine tuning of the magnetic field evolution of the PWN as well as a non-standard 
spin-down history for the pulsar, and it is not clear whether this can be consistently 
applied to other PWN. Alternative models for the acceleration of the radio emitting 
particles at sites in the nebula other than the termination shock should be explored. 

In the above discussion, it has been implicitly assumed that the electrons have 
been injected evenly throughout the nebula, the so-called one-zone model (Clicvalie^ 
200(1) . In a real PWN, the electrons are injected at the termination shock within 
a few tenths of a parsec from the pulsar and then propagate outwards, initially at 
the typical post-shock velocity v ~ c/3 but then with the bulk velocity typically 
decreasing as 1 /r'^ (this is true whatever the emission geometry, since the nebula will 
have an expansion as long as the external confinement is spherically symmetric; 
this is probably not the case for some rP WN - see ti7.6.7|l . If the PWN bubble is large 
enough, the X-ray emitting particles will have a chance to cool significantly before 
reaching the outer edge, and it is expected that the X-ray spectral index will increase 
as a function of distance from the pulsar. For all energies where the travel time to the 
outer edge of the PWN is more than the synchrotron cooling time, then the PWN 
should decrease in size with increasing e nergy. Therefore, we should expect the radio 
PWN to be larger than the X-ray PWN. lKennel fc Coroiiiil l|l984bl) modelled this in 
detail, and found the integrated X-ray spectral index should be Aa = (4+a/9) larger 
than that of the uncooled particles emitting just downstream of the termination 
shock. For reasonable values of a = 0.0 — 1.0, Aa — 0.44 — 0.56, very close to the 
simple, one-zone model estimate of Aa = 0.5. The spatial resolution of Chandra and 
XMM-Newton has all owed this increase in spectral index with radius to be observe d 
in several PWNe Ce.g. lWeisskopf et allEooalSlane etallboOOHPorauet et al.Ll2003|) . 



7.6.5 The Crab and Other Classic "Plerions" 

Until recently, observational and theoretical studies of PWNe have been 
dominated by the Crab nebula. It contains the most energetic pulsar known in the 
Galaxy, it is bright at all wavelengths, has a precisely known age due to historical 
records of the associated supernova, and has been intensely studied for over 100 
years. It had been the defining member of the class of "supernova remnants" with 
no observable shells, sometimes called plerions meaning filled center iWeller, 1978 ), 
and provided the conclusive proof that neutron stars are born in Type II supernovae. 
Ironically, although it has been known as the remnant of SN 1054 A.D. for many 
years, it is something of a misnomer to call it a SNR since the nebula is a result 
of the pulsar wind and not the blast wave of the supernova itself. Why there is no 
observable radio shell is still uncertain, with explanations including that it is in a low 
density region of the ISM and that SN 1054 was an anomalously low-energy event 
I Frail et al. i ll 9951) . Ther e are several other PWNe in the Galaxy with no observable 
shell, typified by 3C 58 llR.evnolds fc Allej. Il985l) . They tend to have much lower 
X-ray luminosites and harder spectra than the Crab. The inferred break energies 
between the radio and X-ray sp ectra of these other P WNe also tend to be much lower 
than that of the Crab's breaks i Woltier et al. . 1997j) . Whether these differences can 
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Fig. 7.14. A composite image of the Crab Nebula showing X-ray in blue, optical 
in green, and radio in red. (Credits: X-ray: NASA/CXC/ASU/J. Hester et ah; 
Optical: NASA/HST/ASU/J. Hester et al.; Radio: VLA/NRAO) (from Chandra 
website http://chandra.harvard.edu) 



be ascribed completely to the difference in spin-down energies is debatable, but 
since observations of the Crab tend to be of much higher quality than of all other 
PWNe, attempting to interpret what is viewed by scaling from the Crab is somewhat 
inevitable. The PWN whose properties are most similar to the Crab's is that around 
PSR B0540— 69, often referred to as the Crab's twin, in the Large Magellenic Cloud 



2001 



I Chanan et all Il984t [Manchester et"al]. Il993t ICotthelf fc Wane 

At radio energies, the Crab appears as an ovoid blob with a superimposed filamen- 
tary structure. The radio spectral index is remarkably spatially uniform ( Bi etcnholz ^^ alJ 



) ^ ^iid the emission is mostly constant excep t for a slow expansion JVehisamv et al 
1992() and a long term overall decline in its flux JAller & R,evnoldsHl985b|r Recently, 




variations in the form of small ripples in the nebula have been detected whose gen- 
eral mor phology correspond to the regions of activity seen in the optical and X-ray 



nebulae ijBietenholz et al.l l200l[) . The radio emission is highly linearly polarized and 



well organized, showing that th e emission mechanism is almost certainly synchrotron 
l|RietRnho1z KronbergilT99nl>. 

At optical energies the Crab is somewhat smaller than in the radio, and consists 
of thermal filaments superimposed on a polarized non-thermal continuum. Near the 
pulsar, there is an underluminous region probably associated with the unshocked 
wind. Outside of this zone, there is a series of bright, non-thermal enhancements, re- 



50 Isolated Neutron Stars 



ferred to as wisps, that have been known for a long time to be variable on timescales 
of months l|Lamplandl. Il92lt lOort fc WalravenL Il956t IScardet Il969|) . Monitoring 
observations with the Hubble Space Telescope have shown the enhancements, the 
innermost one starting at the edge of the underluminous zone which is presum- 
ably at the w ind termination shock, move outwards over several weeks to months 
I Hester et all 11^96) . One model of the radial placement of the wisps posits that 
if the wind contains a significant fraction of ions, the post-shock ions would cause 
enhancements in the magnetic field at multiple s of their c yclotron radius, and hence 
be sites of enhanced synchrotron emission ( Gallant fc Arons. 1994) • 

In X-rays, the Crab has a very definite toroidal morphology with polar jets. Again, 
surrounding the pulsar is an underluminous zone ending in a bright ring with knotty 
enhancements which have been observed with Chand ra to form and move outwards at 
V 0.5c and then diffuse, similar to the optical wisps IjHester 



Outside of 

this region is a thick torus, which ends rather abruptly well within the boundaries of 
the optical and radio nebula (Fig. l7.lHl . There is clear evidence of spectral steepening 
with increasing radius in the torus. The jets emit a much smaller fraction of the X- 
ray luminosity than does the torus, and overall have a harder spectrum. Outward 
motion of features within the jets has also been observed. The proper motion of the 
pulsar is alon g the line of the jets, sugg esting the natal ki ck wa s along the spin axis 
of the pulsar l|Caraveo k Mignanil ^9991 IWeisskopf et allbOOOl) . 

The unpulsed, and hence presumably nebular, emission spectrum is seen to extend 
smoothly from the X-rays all the way to ~ 25 MeV. Above this is a sharp, possibly 

Il99fia|) indicating some limiting factor in the particle 
100 MeV, a rising unpul sed componen t is seen 



Il1. l2002|) . 



variable cut-off I 

acceleration mechanism. Above iuu iviev, a rismg unpui sea coni pc 
which can be observed all the way into the TeV energy range I Weeke£lll99l|) . This is 
most likely due to inverse Compton scattering off the synchrotron emitting electrons 
in the nebula. 

Since the Crab is energetically unique among Galactic pulsars, it is valid to wonder 
if its properties can be legitimately scaled down over two orders of magnitude to other 
pulsars with more typical spin-down energies. 3C 58 was the second plerion to be 
identified and is thought to be a result of SN 1181. Its pulsar, which was only 



recen tly discovered throug h X-ray and very deep radio observations (M urray et al 



l2002l ICamilo et allf2002cl) . has a spin-down energy of 2.7 x 10'^''erg s~^, second only 
to the Crab among Galactic PWNe. Although its spin-down energy is only a factor 
of ~ 15 less than the Crab's and the nebular radio luminosity is a factor of ~ 10 
less than the Crab nebula, the nebular X -ray luminosity is a factor of ^ 1000 less 
( Becker et al .'. '1982: l ^irray et all l2002l) . There appears to be an unusually sharp 
spectral break (,, Green fc Scheueiill992l) at ~ 50 GHz. The radio luminosity appears 
to be incre asing, suggesting the unse en SNR reverse shock may already be affecting 
the PWN jAller fc Revno]d4ll985a|^ . 

The Vela pulsar is generally considered to be the prototype of pulsars with Tc ~ 
10,000 yr and spin-down energies in the 10^^ - 10^"^ erg s ^ range. Its PWN ap- 
pears very bright and can be well-resolv ed at X-ray en e rgies b ecause of its proximity 



to Earth, being only ~ 290 pc away i Dodson et al. . 2003a(l . There is an associ 



ated radio SNR shell, but i t is p atchy, uneven, and diffuse, as is typical of older 
remnants I Clark fc Caswell 19761) . having three sections called Vela X,Y, and Z. 
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Fig. 7.15. 5 GHz ATCA image of the Vela radio lobes near the pulsar with polar- 
ization vectors showing how highly structured the magnetic field is. The contours 
near the center are of the X-ray nebula as seen with Chandra. No te this is only 
a small part of the Vela X radio region. Adapted from Fig. 10 of iDodson et alJ 
J2003bfl 



Weiler fc Panagial l|l98fll) noted that the Vela X region had radio properties more 
similar to the Crab and other plerions than to ordinary SNR shells. Due to the 
age of the Vela SNR, it is likely that the reverse shock encountered the PWN long 
ago, although there may still be some transient effects. The morphology of the Vela 
PWN in X-rays is a double arc with narrow, jet-like enhancements, again near the 
line of proper motion. While this has sometimes been interpreted as a bow-shock 
morphology, it is now usually interpreted as a double torus LHelfand et al.. 2001. 



An_apparently variable, extended jet feature is sometimes prominent ijPavlov et al 
I2OO.1I) outside of the arcs. There is no clear indication of an underluminous region 
near the pulsar, although that may be due to an unfortunate viewing angle, and the 
observed double-arc system may be the location of the termination shock. 

The physical size of the bright part of the Vela X-ray PWN is about an order of 
magnitude smaller than the Crab, and its X-ray efficiency in the 1-10 keV band is 
more than two orders of magnitude less. The X-ray spectrum is flat ter and the high- 
energy cut-off is probably less than 10 MeV I de Jager et al. , 1996b|) . Whether these 
are indicative of a general trend of lower efficiency, spectral hardening, and lower 
maximum particle energies with lower E is diffic ult to say, alth ough there is some 

and the maximum 
EJ-^/^ so the re is a theoretical 
19964 - The X-ray PWN 



evidence for a dependency of spectral index on E I Gottheli Iji 
potential drop across the open field lines Err,.n.:r. 



expectation for a lower high-energy cutoff l|de Jager et al 



is centere d between two radio lobes which are much larger than the X-ray PWN 
CFig. lT.lSfPodson et al. . 2003b(l . and are themselves just a small part near the edge 
of the filamentary, flat spectrum Vela X radio complex. It is likely the reverse shock 
has displaced the larger radio nebula which is why th e proper motion vect or of the 



pulsar is not directed away from the center of Vela X ijBlondin et al.l l200l[l 
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7.6.6 Young Composite SNR 

The standard theoretical picture of a PWN growing within the expanding 
shell of the supernova blast wave has been around since th e early eighties. Althoug h 
quite a few young composite SNR were known at that time ( Helfand &: Beckeil 1987|) . 
in embarrassingly few cases had the exciting pulsar been observed. While the un- 
known radio beaming could be invoked to account for the lack of observable radio 
pulsations, this situation made it difhcult to confront theoretical expectations with 
observational realities. Since the mid-nineties, s everal pulsars have been di scovered 
within young composites by deep radio searches llCamilo et alll200nH2002b^ or with 



sensitive X-ray searches I Torii et al. . 1997t Gotthelf et alll200(]() . representing differ- 
ent stages and scenarios in the early evolution of PWN systems. 

Perhaps the system whose observational properties best match the assumptions 
made in the models is Gil. 2— 0.3. This bright, remarkably circular remnant is often 
assoc iated with a "Guest Star" observed by the Chinese in 386 A.D. l|Clark fc Steohensonl . 

a nd the age inferred from expansion measurements of shell support this asso- 
ciation l|Tam fc Robert4l2003l). HI absorption measurements place the remnant at 
a distance of ~ 5 kpc l|Green et al 111988). A 65 -ms X-ray pulsar is within a few 
arcseconds of the geometrical center of the shell ( Kaspi et al. , 200HJ) which has a 



characteristic age of 24,000 yr; over ten times that of the SNR ijTorii et aL , 1997|) 
strongly implying Pq P and E ~ Eq. High-resolution radio and X-r ay observa- 



tions have separated the P WN from the shell and measured its spectrum l|Tam et al, 
200ilRoberts et ahlEnfii. 



We therefore have in Gil. 2— 0.3 an example of a composite remnant where the 
assumptions of near spherical symmetry of the shell and negligible displacement of 
the pulsar from its birthsite are demonstrably true, the energy output history is well 
constrained to be nearly constant, the age is known, and the distance is fairly well 
determined. Given its age and estimates of the swept-up mass, it should be nearing 
the Sedov phase, so the reverse shock should be nearing the PWN. However, the 
radio PWN is one of t he largest relative to the shell and there is marginal evidence it 
is expanding (|Tam & r!,obertsi . i200.1i) . suggesting the reverse shock has not yet begun 
to crush the PWN bubble. It therefore should still be expanding supersonically, 
and may be shock-heating the surrounding ejecta. There is thermal X-ray emission 
which seems to be morphologically related to the radio PWN, but may also be a 
region of enhanced shell emission seen in projection (Fig. I7.l5|l . The velocity of the 
shell inferred from the thermal emission is consistent with that inferred from the 
expansion rate and HI distance estimates ( Tam fc Robertsl boOSfl . 

Current estimates of the X-ray and radio spectral indices are consistent with a 
single cooling break in between. However, given the age of the nebula, the implied 
magnetic field is unusually high. If it uniformly fills a region of the apparent bubble 
size, then the required magnetic energy is much larger than could have been sup- 
plied by the pulsar. This may indicate that the pulsar wind is highly non-spherical, 
predominantly coming from within ^ 10° of either the poles or the equator. T he 



narrowness of the X-ray emission may support this interpretation IjRoberts et al 



2003). Varying X-ray spots are seen in this nebula as well, with apparent motions of 



much faster than the expected c/3 of a standard post-shock flow. So even with 
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Fig. 7.16. Chandra and VLA image of Gil. 2— 0.3 and its central PWN showing the 
relationship between the soft thermal X-ray emission (red), the radio synchrotron 
emission (green) and the hard, non-thermal X-ray emission (blue). The point 
source at the center is the pulsar, only seen in X-rays. Adapted from Fig. 1 of 
[Roberts et al.. (,2003.1 . 



this most "simple" of PWN/SNR systems where our knowledge is most complete, 
there are anomalies which are not well understood. 

Kes 75 is a young (r < 2000 yr) distant [d - 19 kpc iBecker fc HelfandL Il984 
composite SNR in many ways similar to Gil. 2— 0.3. It co ntains a 324 ms X-ray pulsar 
with the youngest known characteristic age, ^ 700 yr ( Clotthelf et all |2000) . The 
Chandra image of Kes 75 l|Helfand et alJ . l2003bl) of both the sheU and PWN in both 
X-rays and radio of Kes 75 look very similar to Gil. 2— 0.3, and the pulsars powering 
both PWN are very X-ray efficient and have unusually broad profiles. However, on 
closer examination, there are remarkable differences. The X-ray efficiency in the 
0.5-10 keV band of the PWN in Kes 75 is ~ 6.5%, which is the highest known, 
while that of Gil. 2— 0.3 is only ^ 0.2%, a much more typical value. The pulsar's 
spin parameters imply both an unusually large magnetic field (i? ^ 5 x lO^'^ G) and 
that Pq << P, so that the energy input into the PWN has changed dramatically 
over its short history. The implied size of the shell of Kes 75 requires an enormous 
expansion velocity (u ~ 13, 000 km s~ ) and therefore implies a remarkably energetic 
explosion or else the neutron star was born spinning very rapidly and a significant 
fractio n of its rotational kine tic energy was somehow transferred to the expanding 
ejecta l|Helfand et alJ.l2003b|) . 

SNR G292.0-I-1.8 is an example of a slightly older (r ^ 3000 yr) composite system. 
Chandra imaging shows a remnant filled with ejecta mate rial, but ther e is al so a clear 
X-ray PWN slightly offset from the geometrical center 1 Hughes^ et al.ll200lj) . In this 
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Fig. 7.17. Ha image of Black Wido w pulsar PSR B 195 7+20 with Chandra X-ray 
contours. Adapted from Fig. 2 of iStappers et al.l (|200.'j^ . 



case there is no clear morphological separation between the radio PWN and shell 
and it see ms likely that the SNR rev erse shock has recently begun to interact with 
the PWN llGaensler fc Wallacell200.ll) . 



7.6.7 The Duck and Other PWNe with Trail Morphologies 

Several PWNe show a predominantly trail morphology at X-ray, and some- 
times radio, wavelengths extending from an X-ray point source, presumably the 
pulsar, back towards the apparent birth-site. This has been observed in pu lsars both 
within the SNR sheU, such as PSR B1853+01 in W44 jPetre et al.Ll2002l). and out- 
side su ch as coming from PSR B1757— 24, the head of the Duck svstem (|KasDi et aT 
2r)r)lajl . X-ray trails have also been seen from sources with Ha bow-shock neb- 



ulae, such as the Bla ck Widow binary millisecond pulsar system PSR B1957-I-20 



I Stappers et al. . 20031 see Fig. I7.17|l . The apparent radio bow-shock PWNe source 
associated with the possibly variable 7-ray source Gc V J1809 - -2328 also shows a 
clear point source and trail morphology in X-rays ifBraie et al.lboo'^ . These type 
of PWNe are the most common around pulsars with spin-down energies log i? < 36.5 
(in erg s""'^). These nebulae are sometimes referred to as bow-shock nebulae or ram- 
pressure confined nebulae, however it is not clear if either of these terms accurately 
reflect the physical situation in respect to the X-ray and radio nebulae. Therefore, 
we will simply refer to these nebulae as rapidly- moving PWNe (rPWN), where the 
pulsar's motion is rapid relative to the radial expansion of the nebula and probably 
relative to the local sound speed. 

The Duck is an illustrative case. Upper lirn its on the proper motion of the pulsar 
( Gaensler fc Frail 2fl0f]l: iThorsett et al.L l2002l) rule out the X-ray emission being an 
actual trail, since the synchrotron lifetime of the X-ray emitting particles is too short 
for th e pulsar to have mo ved the length of the X-ray emitting region in the required 



time i Kaspi et al. , 2001a(l . The X-ray emission is therefore due to a collimated out- 



flow, somehow constrained to move within a cavity created by the pulsar's passage. 
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A spectacular Chandra X-ray image of the Mouse shows the structure of the X-ray 
nebula to consist of a narrow tail region in a bro ader, fainter, bo w-shaped nebula 
which may be too faint to see clearl y in the Duck llGaensler et all ioOSc). The ex- 
tremely long radio tail of the Mouse ( Predehl fc KulkarnillQQsIl is another example 
where it seems something must continue to confine the flow well beyond where ram- 
pressure from the pulsar's passage sho uld no longer be reley ant. The extraordinary 
jet/trail coming from PSR B1509— 58 l|Gaensler et al. . 2002(l also demonstrates that 
in many cases pulsar outflows can be narrowly confined over large distances even 
while still within the parent SNR. 

The Black Widow pulsar PSR 1957-1-20 l|Fruchter et allllQSSj) . so called because 
its wind is believed to be destroying its white dwarf companion, is interesting for 
several reasons. Unlike the other sources with X-ray trails, it is a millisecond pulsar, 
demonstrating that old, recycled pulsar s also have relativistic winds. It also has 
an Ha nebula ( Kulkarni fc Hesteil Il98 8'l which clearly deline ates the forward bow- 
shock as well as a narrow trail of non-thermal X-ray emis sion dStappers et all l200.l() 
which can serve as a test bed for simulations of rPWNe l|Bucciantinl |2002|) . 



7.6.8 The ^-Ray Connection 

The magnetospherc of a young pulsar is one of the few places in the universe 
where particles can be accelerated such that they will emit high-energy 7-ray emission 
through either the synchrotron process, curvature radiation, or inverse Compton 
scattering. Approximately 1/3 of the known PWNe are coincident with sources 
of emission at £^ > 100 M eV observed by the E GRET instrument on the Compton 
Gamma Ray Ohservatory (Hartm an et al. . 19991) . Some of these are the known 7-ray 
pulsars ( Thompson et al.l.l999j , and this may be an indication that the potential 
needed to produce pulsed 7-ray emission observable at typical Galactic distances 
is similar to that needed to produce observable PWNe. However, several of the 
PWNe are asso ciated with unidentified EGRET sources which appear to be variable 
( Roberts et all EoOl: Nolan et al., 200^ on timescales of a few months, similar to the 
synchrotron cooling timescale at these energies, indicating that the PWN may be the 
source of emission. Practically speaking, 7-ray sources can serve as tracers of PWNe. 
Six of the PWNe listed in Tables 1731 and 17^ were initially discovered through X-ray 
imaging of a 7-ray error box. In addition, the shape of PWNe can be used to infer 
the viewing angle, and possibly the magnetic inclination angle, of the pulsar (e.g. 
iNg fc Roma ni. 2004} which is critical information for testing models of pulsed 7— ray 
emission (^ H7.2(l . Observations from the upcoming GeV energy telescopes ACILE and 
GL AST wiW be able to distinguish between pulsed and unpulsed emission from many 
of these sources, verify any variability, and locate more PWNe buried deep in the 
Galactic plane. 

Steady, unpulsed emission up to several TeV ha s been observed towards three 



pulsars (Crab, Ve la, and PSR B1706-44 IWeekesL Il99lt lYoshikoshi et all Il997t 



Kifune et al. . 199.ll) 1 with PWNe, presumably caused by inverse Compton scattering 
by electrons in the inner nebula. This is direct evidence of very high energy accel- 
eration and puts constraints on the magnetic field in the inner nebula. As the new 
generation of ground based Cerenkov air shower telescopes comes on line, many more 
PWNe will undoubtably be observed at energies of ^ 100 GeV to a few TeV. 
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7.6.9 Current Trends and Future Directions 

The high-resolution X-ray imaging capabihties of the Chandra and XMM- 
Newton telescopes are allowing detailed study of the structure of many PWNe. The 
presence of coUimated jet-like outflows as a common feature of PWNe is a surprise 
which most models of PWNe, dominated by a toroidal equatorial flow, did not allow 
for. New analytical and numerical models are now bei ng d eveloped with v arying 
amount of power being lost through polar outflows JlChangoulian fc Bogovalovt 12003 



Komissarov fc Lvubarskvl 20031: IShibata et all 2003h . PWNe have also been shown 



to be dyna mical X-ray system s with observable variations on timescales as short as a 
few weeks ( Hester et al. , 20021) . As quantitative measurements of these variations are 



made in more PWNe the need for the further development of models of acceleration 
in the inner nebula becomes more acute. 

The remarkable diversity of PWNe has also been a surprising result of the new 
X-ray images. Many of the observations currently archived are too short to distin- 
guish different spectral characteristics of the various features within PWNe. There 
is a suggestive correla tion between th e X-ray spectral index of PWNe and the mag- 
netospheric potential jGotthelAliool . but whether this will be born out by further 



observations and whether there is a different correlation for jet-like and toroidal 
structures will require deep observations of many more PWNe. The structure of 
the PWNe also seems to depend on the magnitude and direction of the pulsar's 
velocity, but so far only a few young pulsars have had their proper motions mea- 
sured. Radio interferometry campaigns now underway (e.g. .Brisken et al.. .2DiKij) 
promise to greatly increase the number of measured pulsar transverse velocities, 
allowing the d etermination of how well jets are aligned with velocity vectors (e.g. 
Dodson et al.l l2003al . if at all. New 2-D and 3-D numerical simulations of PWN 
evolution are beginning to yield intriguing results. For example: relaxing the spher- 
ical symmetry outflow constraint but keeping axial symm etry leads to a multiplicity 
of sh ock regions that can mimic the Crab morphology i Komissarov fc Lvubarskvl 



20031) . Toroidal magnetic outflows can cause an elo ngation of PWNe with spherical 



particle outflows, reminiscent of the shape of 3C58 I van der Swaluwi 20031) . and i 



m- 



teractions with the reverse shock c an lead to c omplicated, asymmetric morphologies 
I Blondin et 3ucciantini et allEoO^ . 



7.7 X-rays from Rotation-Powered Pulsars in Binary Systems 

For completeness, we mention briefly another related source of X-rays from 
rotation-powered pulsars - shock emission in binary systems. Although these systems 
do not obviously satisfy the mandate of this chapter on "isolated neutron stars," 
the nature of their X-ray emission is very similar to that in a PWN. For high-i? 
pulsars in binary systems, the confining medium can be the wind of the companion 
star. In this case, the confinement, particularly for eccentric binaries, is strongly 
orbital-phase dependent, leading to a PWN of luminosity (and, presumably, size) that 
varies systematically with orbital phase. Even circular orbits are of interest since the 
orientation of the PWN, which, in this case, is not expected to be isotropic, changes 
with orbital phase. In principle, such systems are an excellent diagnostic of shock 
acceleration and the pulsar wind as they have regular, repeating and predictable 
dynamical properties. They also represent a new form of "X-ray binary," one not 
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powered by accretion. However, only two such systems have been detected in X-rays: 
PSRs 1259—63 and B1957+20. Their rarity is because of the need to have a high 
E pulsar in a binary system in which the two components come sufficiently close to 
each other for the shock to be strong, but not too close so that accretion occurs. 

The 48-ms radio pulsar PSR B1259— 63 is in a 3.4 yr binary orbit having eccen- 
tricity 0.87 ( Johnston ct al., 1992). From radio timing, the pulsar is known to have 
i? = 3 X 10^^ G and Tc — 3x 10'"' yr. The pulsar's companion, identified by its location 
within the < 1" pulsar timing error box, is the 10th mag B2Ve star SS 2883. It has 
mass ~ 10 Mq and radius ~ 6 R©, deduced from its spectral type. The system pro- 
vides an evolutionary link between the rotation-powered pulsars and the high-mass 
X-ray binaries. For PSR B1259— 63 . near periastron , the p ulsar approaches its com- 
panion to within ~25 Be star radii. Cominskv et alJ 1 1994h detected variable X-rays 
from the system near apastron. This ruled out standard wind accretion scenar- 
ios, m agnetospheric emission, and emission from the companion. ICominskv et al 
l|l994j) suggested either some form of non-standard accretion or the pulsar wind 
shocked by the com panion wind as possible mechanisms for the X-ray emission. 



Tavani et al.l l|1994|) considered pulsar/Be star wind interactions in detail, in par- 
ticular, shock emission at the location of pressure balance between the pulsar and 
Be star winds. They suggested this mechanism produced the apastron X-rays, and 
predicted PSR 1259—63 would be a moderately strong, unpulsed X-ray source near 
periastron, unable to accrete due to the sho ck distance fro m the pulsar being m uch 
larger than the accretion radius. Kaspi et al. (, 1995 ) and Hi ravama et al.l lll99(Jl re- 
ported on ASCA observations around the periastron of 1994; Hiravama et alJ l|l999l) 
summarize all the X-ray observations including two more made at the subsequent 
apastron. The source was clearly seen to increase in intensity by a factor of > 10 
near periastron relative to apastron, with the peak luminosity ^ 10'^^ erg s~^. The 
emission was well described at each epoch by a power law, however the photon index 
clearly varied from 1.6 to 2.0, with the emission softest at periastron. Tavani fc Arond 
( 19971) considered the X-ray radiation mechanisms and interaction geometry in de- 
tail, showing that a synchrotron/inverse-Compton scattering model of emission of 
electron/positron pairs accelerated at the inner shock front of the pulsar cavity and 
adiabatically expanding in the MHD flow explains well the observed time-variable 
X-ray flux and spectrum. They conclude that most likely the Be-star spin axis is mis- 
aligned with the orbital angular momentu m, and its mass o utflow rate was constant 
over the '-~^2-yr period in question. Campana et alJ ()l995f) argue that an unusual 
Galactic X-ray source, LSI-|-61°303, is a similar system to PSR B1259— 63, but with 
the radio pulsations permanently eclipsed. 

PSR B1957-I-20 is a 1.6-ms recycled radio pulsar in a 9-hr binary system with 
a low-mass companion ( Fruchter et al. . 19881) . The pulsar has E — 10'^^ erg s" 



,-1 



For ~10% of every orbit, the radio pulsations are eclipsed by the wind from the 
companion, which is being ablated and, eventually, evaporated, by the pulsar wind. 
For this reason, the pulsar is sometimes known as the "Black Widow," as it seems 
to be destroying the star that gave it new life as a recycled pulsar. The system 
is surrounded by an Ha bow shock nebula, oriented with apex in the direction of 
the pulsar's known proper motion, a result of the interaction o f the pulsar wind 
with the ambient interstellar medium ( Kulkarni fc Hest'eR 1988j) . Arons fc Tavanil 
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l|1994|) argued that intra-binary shock emission ought to be observable, because of 
the interaction of the pulsar wind with material being ablated off the companion. 
They also argued that X-rays might originate from the bow shock. Of course, mag- 
netospheric X-rays might be detectable from this pulsar as well. The first X-ray 
detection of PSR B1957+20 came from R OSAT, althoug h only handful of photons 
were detected (^ruchtcr ct al., 1992; Kulka rni et al. . ^"^^ the conclusions that 



could be drawn were limited. More recently, Chandra observed the PSR B1957+20 
system (see Fig. I7.17II . and detected both a point source as well as an extended X- 
ray tail (.Stappers et a l.. 2003). The tail emission is likely related to the bow shock 
nebula, and represents the first proof that millisecond pulsars have relativistic winds 
like their much younger counterparts. The origin of the point-source emission is 
less clear as the observation's time resolution was insufficient to detect pulsations. 
Planned XMM-Newton observations should be able to decide with certainty whether 
the point source is due to magnetospheric emission, or if it is from the intra-binary 
shock, as in PSR B1259-63. 

The authors thank B. Gaensler and D. Yakovlev for discussions, for careful reading 
of the manuscript and for many helpful comments. We also thank P. Jaikumar, D. 
Lai, M. Lyutikov, G. Pavlov, S. Ransom, R. TuroUa, M. van Kerkwijk, and S. Zane 
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